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We  present  results  from  an  experimental  study  of  the  decoherence  and  decay  of  quantum  states  of  a  trapped 
atomic  ion’s  harmonic  motion  interacting  with  several  types  of  engineered  reservoirs.  We  experimentally 
simulate  three  types  of  reservoirs:  a  high-temperature  amplitude  reservoir,  a  zero-temperature  amplitude  res¬ 
ervoir,  and  a  high-temperature  phase  reservoir.  Interaction  with  these  environments  causes  the  ion’s  motional 
state  to  decay  or  heat,  and  in  the  case  of  superposition  states,  to  lose  coherence.  We  report  measurements  of  the 
decoherence  of  superpositions  of  coherent  states  and  two-Fock-state  superpositions  into  these  reservoirs,  as 
well  as  the  decay  and  heating  of  Fock  states.  We  confirm  the  theoretically  well-known  scaling  laws  that  predict 
that  the  decoherence  rate  of  superposition  states  scales  with  the  square  of  the  “size”  of  the  state. 

PACS  number(s):  42.50.-p,  03.67.Lx,  03.65.-w,  32.80.Pj 


I.  INTRODUCTION 

The  decoherence  and  decay  of  quantum  states  coupled  to 
a  variety  of  reservoirs  have  been  investigated  extensively  in 
theory  [1-8].  The  model  in  these  studies  is  typically  a  sys¬ 
tem  harmonic  oscillator  coupled  to  a  bath  of  environment 
harmonic  oscillators.  One  of  the  most  interesting  results  has 
been  the  realization  that  macroscopic  superposition  states  de¬ 
cay  at  extremely  fast  rates.  As  an  illustration,  consider  a 
charged  harmonic-oscillator  system  in  a  superposition  of  co¬ 
herent  states  i/f«|  cr [ )  + 1  a2)  coupled  to  the  noisy  electric 
field  E  of  the  environment.  The  interaction  potential  is  V 
=  —  qx  ■  E,  where  x  is  the  position  and  q  is  the  charge  of  the 
particle.  Regardless  of  the  temperature  of  the  reservoir, 
(E2)  A0,  so  that  the  system  is  always  subject  to  some  level 
of  noise  from  the  environment.  For  a  superposition  of  coher¬ 
ent  states  coupled  to  a  reservoir  of  fluctuating  fields,  a  simple 
scaling  law  may  be  stated  [3]:  the  rate  of  decoherence  scales 
as  the  square  of  the  separation  of  the  wave  packets,  A  aj 2 
=  \al  —  a2\2.  In  an  idealized  case,  the  superposition  is  cre¬ 
ated,  the  amplitude  reservoir  is  coupled  to  the  system  for  a 
time  f,  and  the  coupling  is  then  turned  off.  The  remaining 
coherence  C  between  the  two  wave  packets,  expressed  for 
instance  as  the  magnitude  of  the  off-diagonal  components  of 
the  system  density  matrix,  is 


C(f)  =  C(0)exp 


A  a 


(1) 


where  y  is  a  coupling  constant  between  the  reservoir  and  the 
system.  For  the  same  interaction,  the  energy  H  of  the  system 
decays  like 


H(t)  =  H(0)e\p[  —  yt],  (2) 
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at  a  rate  independent  of  the  size  of  the  initial  state.  For  su¬ 
perpositions  of  macroscopic-sized  wave  packets,  the  quan¬ 
tum  coherence  may  be  obliterated  in  a  time  over  which  the 
energy  of  the  system  does  not  change  appreciably. 

Significant  attention  has  been  paid  to  the  role  of  decoher¬ 
ence  in  the  classical/quantum  correspondence  and  in  funda¬ 
mental  issues  of  physics.  Recently,  decoherence  has  also 
been  studied  in  a  more  pragmatic  role:  as  a  primary  impedi¬ 
ment  to  quantum  computing.  Quantum  computation  relies  on 
entanglement  of  large  quantum  states  to  perform  efficient 
calculations.  Such  states,  separated  by  large  distances  in  Hil¬ 
bert  space,  will  be  very  susceptible  to  the  detrimental  effects 
of  decoherence.  Any  system  will  interact  with  its  environ¬ 
ment,  so  the  nature  of  the  interaction  and  the  time  scales  over 
which  it  acts  are  of  critical  import  in  implementations  of 
quantum  information  processing.  Since  trapped  ions  are  a 
leading  technology  for  investigations  in  quantum  computing, 
a  study  of  decoherence  in  an  ion  trap  quantum  computer 
system  [9],  as  presented  here,  is  particularly  relevant. 

Over  the  past  several  years,  techniques  have  been  de¬ 
scribed  to  generate  mesoscopic  superpositions,  called 
“Schrodinger  cats,”  of  motional  states  of  trapped  ions  [10] 
and  of  photon  states  in  cavity  quantum  electrodynamics  [1 1], 
In  both  cases  decoherence  through  coupling  to  ambient  res¬ 
ervoirs  was  observed,  and  sensitivity  to  the  size  of  the  super¬ 
position  was  demonstrated  [10,1 1].  In  more  recent  work  with 
trapped  ions  [12],  we  extended  the  investigations  beyond 
ambient  reservoirs  and  “engineered”  the  form,  bandwidth, 
and  strength  of  coupling  to  the  reservoir.  We  observed  the 
quantitative  scaling  laws  of  the  decoherence  of  superposi¬ 
tions  of  coherent  states  and  superpositions  of  Fock  states  into 
a  variety  of  reservoirs  for  a  range  of  parameters.  This  paper 
is  a  more  detailed  discussion  of  the  experiments  of  Ref.  [12] 
with  additional  theoretical  development  and  additional  data. 

The  paper  is  arranged  as  follows.  In  Sec.  II  we  review  the 
theory  of  the  damping  of  a  harmonic  oscillator  coupled  to 
several  types  of  reservoirs.  In  particular,  we  derive  specific 
formulas  for  the  time  evolution  of  the  density  matrix  for 
several  types  of  system-reservoir  interactions.  In  Sec.  Ill  we 
discuss  the  experimental  apparatus  and  techniques,  and  in 
Sec.  IV  we  analyze  the  data. 
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H.  THEORETICAL  SUMMARY 
A.  Damping  of  harmonic  oscillators 
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1.  Thermal  amplitude  reservoir 

We  begin  by  reminding  the  reader  of  some  formal  results 
in  the  theory  of  damping  of  harmonic  oscillators  coupled  to 
two  types  of  baths.  In  this  section  we  consider  the  case  in 
which  the  reservoir  is  a  bath  of  oscillators,  and  the  coupling 
is  to  the  position  of  the  oscillator  via  the  interaction  Hamil¬ 
tonian 
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where  pnm( 0)  is  the  initial  density  matrix  and 
N(t)  =  n(  \-e~yt) 


(7) 

(8) 


where  T,  is  the  coupling  rate  of  the  /th  bath  oscillator  to  the 
system,  is  the  lowering  operator  of  this  reservoir  oscilla¬ 
tor,  and  a  is  the  lowering  operator  of  the  system  oscillator. 
This  is  known  as  an  amplitude  coupling.  If  furthermore,  the 
reservoir  oscillators  are  in  thermal  equilibrium  at  tempera¬ 
ture  T,  then  the  system  exchanges  energy  with  the  reservoir 
in  a  process  leading  to  overall  thermal  equilibrium.  The  time 
evolution  of  this  process  is  described  by  the  following  master 
equation  for  the  system  density  matrix  formed  by  tracing 
over  the  reservoir  degrees  of  freedom  [5]: 


y  _  A  A  A 

p(t)=  —  (n  +  l)(2flpfl*—  a^ap  —  pa  *a) 


y — 

+  — n{2a  ^  pa  —  aa  1  p  —  paa  ' ), 


(4) 


where  y  is  the  system  decay  rate  and  n  is  the  average  num¬ 
ber  of  quanta  in  the  reservoir  at  the  resonance  frequency  of 
the  system  oscillator  a>0 : 


can  be  interpreted  as  the  average  number  of  quanta  in  the 
system  at  time  t  for  an  initial  ground-state  system  [p(0) 
=  1 0 ) (0 1 ] .  The  form  of  Eq.  (7)  is  qualitatively  illuminating: 
the  evolution  of  a  given  density-matrix  element  depends  only 
on  neighbors  along  its  own  diagonal.  We  note  that  this  cal¬ 
culation  has  been  performed  in  Refs.  [13]  and  [14].  We  used 
Ref.  [13]  to  derive  Eq.  (7);  in  this  reference  the  final  form 
[Ref.  [13],  Eq.  (37)]  seems  to  be  written  incorrectly,  though 
the  correct  result  can  be  derived  from  earlier  equations.  In 
Ref.  [14]  the  result  is  written  in  terms  of  the  Q  function,  and 
agrees  with  Ref.  [13],  and  therefore  not  with  our  result. 

In  some  cases,  we  are  interested  in  the  time  evolution  of 
the  diagonal  elements  only,  which  are  given  by: 

1  ^  /  N(t)  yi  e ~ 7,12  \2n~2j 

Pnn(t)  -  i  +N(t)  jit)  {  1  +  A T(t)  )  \  1  +A T(t)  )  ^ 

"  /  e-y  W n  +  l-j\ 

^o'1  1  +N(t)J  \  n-j  ) 

X(jJPn  +  l-j,n  +  l-j(0)-  (10) 


e~hio0lkBT 
l  —  e~f‘a>0/kBT 


As  stated  above,  the  zero-temperature  reservoir  is  a  spe- 
(5)  cial  case  of  the  amplitude  reservoir.  We  write  the  solution  in 
this  case  as: 


and  kB  is  Boltzmann’s  constant.  The  zero-temperature  reser¬ 
voir  (n  =  0)  is  a  special  case  of  the  thermal  amplitude  reser¬ 
voir  in  which  the  inescapable  quantum  noise  dominates  the 
classical  thermal  noise  of  fluctuating  oscillators. 

a.  Amplitude  damping,  Fock  states.  The  solution  of  Eq. 
(4)  is  nontrivial,  and  has  been  discussed  at  some  length  in  the 
literature  (which  is  reviewed  in  Ref.  [8]).  We  first  consider 
its  solution  in  a  Fock-state  basis.  In  this  basis,  the  decompo¬ 
sition  of  p(f)  is  given  by 

p(0  =  2  Pnm(t)\n)(m\.  (6) 

n,m 

Based  on  calculations  presented  in  Ref.  [13],  we  state  the 
solution  to  Eq.  (4)  with  the  expansion  of  Eq.  (6)  as: 
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Pnn{t)  =  e-n^  (l-e-*)]  pB+u+/(0).  (12) 

1=0  \  n  I 

b.  Amplitude  damping,  coherent  states.  In  this  paragraph 
we  consider  the  solution  of  Eq.  (4)  in  a  situation  appropriate 
to  coherent-state  evolution.  We  follow  exactly  the  results  of 
Ref.  [4].  For  coherent  states  the  density  matrix  is  expanded 
as: 
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p(t)=  J  d 2 lid 2vP(n,v*)  ,  |  '  ,  (13) 

where  P(p,,v*)  is  the  positive-P  representation  [5-7].  The 
solution  to  Eq.  (4)  is  written  [4] 

p{t)=  J  d2fid2vP0(/i,v*) 

X{v\fi)-e~^1+2N^t{J]i)pT(t)t>\v),  (14) 

where  P0(/i,v*)  describes  the  initial  system  density  matrix 
through  Eq.  (13), 

e  ~  y"2 

p=  [+2N(t)m+N(t)]p  +  N(t)v},  (15) 

e-yin 

v=  1+2jV(f>{[1+^(f)]^+^(l)^},  (16) 


and  Pr(t)  is  a  thermal  state  density  matrix  with  N(t)  quanta: 


Pr(t)  = 
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(17) 


N(t )  is  defined  in  Eq.  (8).  V  is  the  displacement  operator 
A«)|o)=|  a).  The  form  of  this  solution  makes  it  difficult  to 
see  the  simple  result  for  the  off-diagonal  elements  of  the 
density  matrix  as  presented  in  the  Introduction.  The  appro¬ 
priate  experimental  measurement  will  reveal  it  when  we  ap¬ 
ply  this  formal  result  below. 

2.  Phase  reservoir 

Here  we  consider  a  reservoir  whose  coupling  preserves 
the  energy  of  the  system,  with  interaction  Hamiltonian 


Hj=h E  (ri£1-a^a  +  r*£ja^a).  (18) 

i 


This  is  known  as  a  phase-damping  reservoir. 

For  the  system-reservoir  interaction  described  by  Eq.  (18) 
and  a  reservoir  in  thermal  equilibrium  at  a  temperature  cor¬ 
responding  to  a  mean  occupation  number  n  [as  in  Eq.  (5)], 
the  master  equation  is  given  by  [6]: 

p{t)  =  —  [2a  ' apa^ a  —  (a  1  a)2 p  —  p(a  'a)2].  (19) 

Here  K  =  K(2n  + 1),  is  the  system  decay  rate  with  K  the 
coupling  to  the  reservoir  of  oscillators. 

a.  Phase  damping,  Fock  states.  For  initial  Fock  states,  the 
solution  of  Eq.  (19)  is  straightforward.  Using  the  expansion 
of  Eq.  (6), 


Pnm=~  2  (n~mYPn 


(20) 


with  solution 
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pnm(t)  =  e-0d2Xn-m)2,pmn{  0).  (21) 

Here  it  is  clear  that  the  off-diagonal  density-matrix  elements 
decay  exponentially  at  a  rate  proportional  to  the  square  of 
their  distance  from  the  diagonal. 

b.  Phase  damping,  coherent  states.  For  the  coherent  state 
basis,  the  solution  of  Eq.  (19)  is  rather  complicated.  We  will 
forego  its  solution  in  favor  of  the  technique  presented  in  the 
next  section. 


B.  Interferometry 

This  section  presents  an  experimentally  motivated  ap¬ 
proach  to  calculating  the  decay  of  quantum  coherences.  Here 
the  off-diagonal  matrix  elements  are  measured  as  the  contrast 
of  an  ion  interferometry  experiment.  This  section  is  adapted 
from  Ref.  [15]. 


1.  Ion  states  and  transitions 

First  we  review  a  few  details  necessary  to  understand  the 
ion  interferometer  analysis.  For  the  purposes  of  this  section, 
a  trapped  ion  has  an  external  degree  of  freedom  that  is  a 
perfect  one -dimensional  (ID)  harmonic  oscillator  (along  z) 
of  frequency  to.  The  internal  degree  of  freedom  is  assumed 
to  be  a  perfect  two-state  system,  the  states  of  which  we  label 
in  analogy  with  a  spin- 1/2  system  by  ||)  and  |]~).  The  states 
1 1 )  and  |t)  are  coupled  to  the  motional  harmonic-oscillator 
states  | n)  via  a  classical  optical  laser  field.  The  interaction 
Hamiltonian  has  resonant  couplings  between  internal  states 
Is)  (ls)=ll)  or  It))  and  motional  states  | n)  whose  matrix 
elements  are  [16,17] 


(s' ,n  '\HI\s,n)  =  h£l(s' ,n '  |  cr+e,,7(fl+at) 

+  o-_e-i’(“+“t)|  s,n),  (22) 

where  cr+  (cr_)  is  the  atomic  raising  (lowering)  operator,  a1 
(a)  is  the  harmonic-oscillator  creation  (annihilation)  opera¬ 
tor,  fl  is  the  laser-atom  coupling  strength,  or  Rabi  frequency, 
and  r]  =  Sk  zo  is  the  Lamb-Dicke  parameter,  where  Sk  is  the 
z  component  of  the  applied  field  wave  vector  and  z0 
=  yJh/2mto  is  the  harmonic-oscillator  characteristic  length 
[16]. 

The  excitation  spectrum  of  the  ion  consists  of  a  carrier 
transition  (||)|n)<->|T)|n))  and  motional  sideband  transi¬ 
tions  (|].)|n)<->|T)|w)n#w).  The  sideband  transition  fre¬ 
quencies  differ  from  that  of  the  carrier  by  a  frequency  of 
(n-m)to,  where  n  —  m  is  the  order  (first,  second,  etc.)  of 
the  sideband.  With  resolved  sideband  cooling  [18]  and  opti¬ 
cal  pumping  we  prepare  the  ion  in  the  ||)|0)  state  before 
each  repetition  of  the  experiment. 

At  the  end  of  any  experiment,  there  is  one  observable:  the 
probability  that  the  ion  is  found  in  the  |  I. )  state,  denoted  by 
P  |  .  The  measurement  is  not  sensitive  to  the  motional  state 
of  the  ion,  so  a  trace  over  the  motion  is  required  to  calculate 
this  probability.  The  final  density  matrix  at  the  end  of  an 
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experiment  is  written  pf,  and  in  general  is  not  separable 
between  spin  and  motion  parts  [19].  The  detection  signal  is 
given  by 

^  =  U|TrMP/||>  =  2  <«Kl|p/|l)|«).  (23) 

n 

For  pure  states  pf=  \  <A/)(  lA/'l>  this  reduces  to 

Pl  =  2  KlK«l<A/)|2-  (24) 

n 

In  experiments  in  which  a  coherence  is  measured,  we  em¬ 
ploy  a  type  of  Ramsey  interferometry.  Genetically,  the  pro¬ 
cedure  is  to  embed  the  coupling  to  the  reservoir  between  the 
initial  and  final  “beam  splitters”  of  a  Ramsey  interferom¬ 
eter.  The  beam  splitters  create  and  undo  a  motional  quantum 
state.  As  the  coherence  of  the  quantum  state  is  degraded  (due 
to  interaction  with  the  reservoir)  the  contrast  of  the  Ramsey 
fringes  decreases,  giving  a  measure  of  the  remaining  coher¬ 
ence. 

2.  Amplitude  reservoirs 

a.  Schrodinger  cat  interferometer  and  amplitude  reser¬ 
voir.  The  first  type  of  interferometer  that  we  employ  uses  the 
Schrodinger  cat  state  [10] 

k)=^(U)K)  +  lT)l«T>),  (25) 

where  ||)  and  |J,)  are  internal  states  of  the  ion,  and  |ajt) 
are  coherent  states  of  the  ion’s  motion  in  the  harmonic  well 
of  the  trap.  As  will  be  seen,  the  action  of  an  amplitude  res¬ 
ervoir  on  this  Schrodinger  cat  interferometer  clearly  illus¬ 
trates  the  scaling  of  the  decoherence  as  a  function  of  the  size 
of  the  superposition. 

A  diagram  of  the  Schrodinger  cat  interferometer  is  pre¬ 
sented  in  Fig.  1.  Initially  the  ion  is  prepared  in  the  state 
|  i/f)  =  1 1)|0).  Next,  we  generate  the  spin  superposition  \f) 
=(ID  +  lt))|oW2,  as  shown  in  panel  1.  The  motional  su¬ 
perposition  is  created  by  applying  a  spin-dependent  optical 
dipole  force  that  approximates  the  coherent  displacement  op¬ 
erators  [10,19]  2?(a|)  and  T>(af)  acting  on  the  states  1 1)  and 
|t).  The  resulting  state  (with  constant  phase  factors  sup¬ 
pressed),  shown  in  panel  2,  is  that  of  Eq.  (25).  These  two 
operations  constitute  the  first  “beam  splitter,”  or  Ramsey 
zone  of  the  interferometer.  During  the  Ramsey  waiting  time, 
the  amplitude  reservoir  [with  coupling  shown  in  Eq.  (3)]  is 
applied,  by  placing  noisy  fields  on  the  trap  electrodes  (this 
will  be  discussed  in  Sec.  Ill  C  1  below).  The  application  of 
the  amplitude  reservoir  for  a  fixed  time  results  in  a  random 
displacement  f3  that  displaces  both  spin  states  equally.  Since 

D(\)|  /<•)  =  e~'Im[x*'c]|A.  +  k),  (26) 

where  |  k)  and  |  A.  +  k)  are  coherent  states,  after  application 
of  the  operator  V{/3)  to  the  wave  function  from  Eq.  (25) 


1)  tt/2  on  spin 

2)  displacement 
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FIG.  1.  Pictorial  representation  of  the  Schrodinger-cat  interfer¬ 
ometer  with  coupling  to  an  amplitude  reservoir.  The  coherent  states 
are  represented  by  minimum  uncertainty  circles  in  phase  space, 
while  the  different  spin  states  are  indicated  by  the  hatching.  In  panel 
1 ,  a  spin  superposition  |i/f)  =  (|j,)-l-|t))|0)/\/2is  created  with  a  ir/2 
pulse  on  the  carrier  transition.  The  cat  state  |</'}  =  (li)lQD 
+ 1|  )|  oq})/\/2  is  completed  with  a  spin  dependent  dipole  force,  as 
shown  in  panel  2.  The  effects  of  coupling  to  an  engineered  ampli¬ 
tude  reservoir  are  shown  in  panel  3.  A  random  displacement  by  f3  in 
panel  3  displaces  both  spin  states  equally,  but  with  different  phase 
factors.  After  the  reversal  of  the  initial  displacement  in  panel  4  the 
motion  may  again  be  factored  out,  but  leaving  a  relative  phase  shift 
2  Im[/3Ao:*]  between  the  spin  components.  An  ensemble  average 
over  y8  leads  to  a  loss  of  contrast  (decoherence)  that  scales  with  the 
size  of  the  superposition  |  A o-| . 

l^)=i(ll)ki  +  t8)  +  enm[M“*1|T)|«T  +  y8»,  (27) 
V  2 

where  Aa=aj-aj  and  an  overall  phase  factor  has  been 
removed.  The  motional  superposition  is  reversed  by  applying 
and  resulting  in 

|«A>=  ^(U)  +  ^2iIm[^A“*]|t»l/S>-  (28) 

A  final  7r/2  pulse  transforms  the  spins  as 

+  (29) 

|TH^=(|T>-eil»,  (30) 

V2 

and  yields  the  state 

\f)=^[(l+ei^s+21m[<SAa*]))\l) 

+  (e2im^a*]-eis)\J[)]\/3),  (31) 

where  S  is  the  Ramsey  phase,  i.e.,  the  phase  difference  be¬ 
tween  the  initial  and  final  it  12  spin  flips.  These  last  two  steps 
constitute  the  second,  recombining  beam  splitter  of  the  inter¬ 
ferometer. 

In  a  single  experiment  the  probability  that  the  ion  is  in  the 
||)  state  is 
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P  |(/3)  =  1(1 1  *A)|2  =  y  C1  +cos(<5+2  H^Ao!*])]. 

(32) 

If  [3  were  fixed  on  each  repetition  of  the  experiment,  as  the 
Ramsey  phase  S  is  swept,  the  signal  would  be  a  perfect  co¬ 
sine  with  phase  dictated  by  the  displacement  f3.  The  deco¬ 
herence  (and  loss  of  contrast  in  this  picture)  enters,  because 
j3  is  a  random  variable  that  fluctuates  to  a  different  value  on 
each  repetition  of  the  experiment,  giving  rise  to  a  random- 
phase  shift  of  the  fringe  that  averages  the  contrast  away  from 
unity.  In  the  experiment,  we  force  (3  to  be  a  Gaussian  dis¬ 
tributed  random  variable,  as,  for  example,  realized  by  a  hot 
resistor  coupled  to  the  trap  electrodes.  We  take  both  Re  (3 
and  Im  j3  to  have  a  standard  deviation  cr.  An  ensemble  aver¬ 
age  of  P  j(/3)  yields  the  form  of  the  Ramsey  fringes: 

P[  =  2  [1  +e~2lA“l~°’2cos  £>]■  (33) 


We  see  that  the  contrast  of  the  fringes  is  exponentially  sen¬ 
sitive  to  (the  square  of)  the  “size”  |Aa|  of  the  superposi¬ 
tion.  This  scaling  [3,7,20]  is  one  of  the  key  results  that  is 
demonstrated  experimentally  in  a  later  section. 

(b)  Fock  state  interferometer  and  amplitude  resen’oir.  A 
second  type  of  interferometer  generates  a  superposition  of 
two  Fock  states  with  the  first  beam  splitter  using  the  tech¬ 
nique  described  in  Ref.  [17].  For  example,  consider  the  in¬ 
terferometer  that  generates  the  superposition  (\0)  +  \2)) /  \f2 
and  then  recombines  the  two  parts  of  the  wave  function  to 
produce  interference  fringes.  Starting  from  the  state  ||)|0),  a 
7t/2  pulse  on  the  first  blue  sideband  generates  the  state  |  f) 
=  (l-l)|0)  +  |T )|  1))/ \/2^-  A  7 t  pulse,  tuned  to  the  first  red 
sideband,  drives  the  transition  |t)|l)— >||)|2)  while  leaving 
the  population  in  the  ||)|0)  state  unperturbed.  This  combi¬ 
nation  of  pulses  generates  the  superposition 

!‘A>=^llXI°>  +  |2».  (34) 

The  second  beam  splitter  of  the  interferometer  is  realized  by 
reversing  the  above  procedure — driving  a  tt  pulse  on  the  first 
red  sideband  followed  by  a  7r/2  pulse  on  the  first  blue  side¬ 
band.  When  the  phase  of  the  second  7r/2  pulse  is  swept,  the 
spin  state  at  the  end  of  the  experiment  oscillates  between  |  [ ) 
and  |t). 

For  the  experiments  reported  here,  we  use  two  pulses  to 
generate  the  beam  splitters  in  the  Fock  state  interferometer. 
For  the  first  beam  splitter,  the  first  pulse  is  a  tt/2  pulse  on  the 
n '  th  order  blue  sideband,  generating  the  state 

l‘A)=^(ll>l°)  +  lT)l»'))-  (35) 

A  second  pulse  selectively  drives  only  one  of  the  spin  com¬ 
ponents,  resulting  in  the  Fock  state  superposition 
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TABLE  I.  The  stimulated  Raman  pulses  required  for  the  first 
beam  splitter  of  the  Fock  state  interferometer  where  the  initial  state 
is  always  ||)|0).  To  reverse  the  superposition  in  the  second  beam 
splitter,  first  the  it  pulse  is  applied,  followed  by  the  tt/2  pulse  with 
a  phase  added  with  respect  to  the  first  ttI2.  RSB  and  BSB  stand  for 
red  sideband  and  blue  sideband,  respectively. 


tt/2  pulse 

77  pulse 

State 

carrier 

lsl  RSB 

||)(|0)  +  |1»/V2 

1st  BSB 

lsl  RSB 

||>(|0)  +  |2))/V2 

2nd  BSB 

lsl  RSB 

U>(|0)  +  |3»/v/2 

1 st  BSB 

2nd  BSB 

|T)(|1>  +  |2»/V2 

I<A> 


l 

n/2 


s)(\m)  +  \n)), 


(36) 


where  |s)  denotes  one  of  the  spin  states  ||)  or  ||).  The  four 
pairs  of  Fock  states  generated  in  the  experiments,  along  with 
the  pulses  used  to  generate  the  superpositions,  are  shown  in 
Table  I. 

The  Fock  state  interferometer  can  be  coupled  to  an  am¬ 
plitude  reservoir.  For  simplicity,  we  restrict  ourselves  to 
Fock  state  superpositions  of  the  form 

l‘A>=^ll>(l°>  +  l,»»-  (37) 


It  is  an  extended  exercise  to  calculate  results  for  the  entire 
interferometer  in  the  manner  of  the  preceding  case  (coherent 
states  in  amplitude  reservoir).  The  complication  is  twofold: 
First,  displacing  Fock  states  leads  to  cumbersome  expres¬ 
sions.  Second,  there  is  a  problem  with  the  beam  splitters:  the 
pulses  of  the  second  beam  splitter  are  tt  and  tt/2  pulses 
strictly  for  the  transitions  1 J, )  |  m )  <->  |  T  )  |  m  —  1 )  and 
ll)|0)  <->  1 1  )\m  —  1 ),  respectively.  Thus  applying  these 
pulses  to  an  ion  in  a  motional  state  that  has  changed  from  its 
original  form  does  not  result  in  full  spatial  overlap.  It  turns 
out  that  while  the  calculation  is  complicated,  the  details  are 
not  important  to  the  actual  experiment,  which  is  performed 
only  for  small  m,  and  in  the  limit  yt  <§  1 .  Rather  than  present 
the  full  calculation,  we  will  merely  state  a  simple  result 
based  on  the  master-equation  approach.  For  states  of  the  type 
shown  in  Eq.  (37),  the  time  evolution  of  the  relevant  density- 
matrix  elements  predicted  by  Eq.  (7)  reduces  to  the  simple 
form 


1 


P0,m 


2(1  +  nyt)1+m 


(38) 


for  yt<  1 .  The  measured  contrast  is  simply  twice  this  matrix 
element.  For  the  range  of  m  and  n  yt  studied  in  the  experi¬ 
ment,  this  expression  differs  by  less  than  5%  from  an  exact 
treatment  of  the  imperfect  spatial  overlap. 
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FIG.  2.  Pictorial  representation  of  the  Schrodinger-cat  interfer¬ 
ometer  with  coupling  to  a  phase  reservoir.  The  coherent  states  are 
represented  by  minimum  uncertainty  circles  in  phase  space,  while 
the  different  spin  states  are  indicated  by  the  hatching.  In  panel  1,  a 
spin  superposition  | <A)  =  (| !)  +  |T))|0)/\/2  is  created  with  a  tt/2 
pulse  on  the  carrier  transition.  The  cat  state  | f)~  (1 1)1  a[) 
+  |T)|<*f))/\/2  is  completed  with  a  spin  dependent  dipole  force,  as 
shown  in  panel  2.  Coupling  to  a  phase  reservoir  is  depicted  in  panel 
3.  The  random-phase  shift  in  panel  3  prevents  the  correct  reversing 
of  the  initial  creation  of  the  cat,  resulting  in  a  loss  of  contrast  due  to 
both  the  phase  shift  f  as  well  as  the  reduced  spatial  overlap. 

3.  Phase  reservoirs 

(a)  Schrodinger  cat  interferometer  and  phase  reservoir.  A 
second  reservoir  that  we  apply  is  a  phase  reservoir.  A  phase 
reservoir  is  realized  by  adiabatically  changing  the  trap  fre¬ 
quency  for  a  period  of  time  and  then  restoring  the  original 
trap  frequency.  If  the  deviation  of  the  trap  frequency  is 
Sco(t),  then  the  phase  shift  of  the  ion  oscillation  in  time  T  is 

rr 

4>=  8a>(t)  dt.  (39) 

Jo 

The  potentials  required  to  change  the  trap  frequency  are  ap¬ 
plied  to  the  trap  electrodes  during  the  Ramsey  time  between 
the  beam  splitter  pulses  (see  Sec.  Ill  C  2). 

The  effect  of  a  phase  reservoir  on  a  Schrodinger  cat  state 
is  diagrammed  in  Fig.  2.  The  generation  of  the  Schrodinger 
cat  state  ]i/f)  =  (|J.)|o'|)  +  |t)|a|))/\/2  is  accomplished  as 
discussed  previously.  A  phase  shift  of  the  ion  oscillation  by 
the  random  variable  <p  introduces  complex  phase  factors  to 
the  coherent  states, 


l<A)=^(ll)l«ie^)+|t)|aTe^)).  (40) 


In  a  single  run  of  the  experiment,  the  probability  to  find  the 
ion  in  the  |  [ )  state  is 

/Ji(^)  =  Kll(A)|2=  y{l-exp[-|Aa|2(l-cos  </>)] 

Xcos[<5+  (|  —  |  a^|2)sin  <J> 

+  2(1  —  cos  <J>)Im  a*a^]}.  (42) 

The  phase  shift  <J>  is  taken  to  be  a  Gaussian  distributed 
random  variable  with  standard  deviation  cr.  The  average  of 

Pf(t 0  is 

p‘4{ 1  “  ~  T7H4"|J(1  “c“'w) 

Xcos[(|  aj2  —  |  a||2)sin  cf>] 

Xcos[<5+2(  1  —  cos  <J>)Im  a*a^]dtf>  L  (43) 


where  the  odd  part  of  the  integrand  vanishes.  We  note  that  in 
general  there  is  a  phase  shift  to  the  fringes  even  though  the 
distribution  of  <fi  is  symmetric.  In  general  this  integral  does 
not  simplify.  However,  in  the  case  of  experimental  interest, 
(specifically,  oq=  —  2af)  in  which  case  Im  a*  af| 
=  0.  We  further  make  the  small-angle  approximation  in  the 
integrand,  which  is  valid  if  cr2<?  1,  when  the  term  e  ll>  /2<r 
dominates;  or,  if  |  A  o'!  1  and  cr^rr,  when  the  term 

g-|Aa|  cos  (f>)  dominates.  In  these  regimes,  the  integral  can 
be  evaluated  analytically  to  yield  the  signal 


II 

cos  8 

-(KI2-H2)V 

1 

2| 

l1  v/l  + 1  Aa]2cr2  CXP 

2(1  +  |Ao'|2cr2) 

1 

(44) 

This  approximate  expression  provides  a  good  guide  to  the 
behavior  and  can  be  used  for  a  simplified  comparison  to  the 
data.  We  note  that  in  the  experimentally  accessible  regime  of 
|Aa|~l-4  and  cr^2,  Eq.  (44)  is  accurate  to  about  10% 
even  though  the  small  angle  approximation  is  not  strictly 
valid. 

(b)  Fock  state  interferometer  and  phase  reservoir.  We 
start  in  the  state  [|«)  +  |m)]|J.)/y/2.  We  again  denote  the  ion 
oscillation  phase  shift  by  cf> ,  and  the  phase  shifted  superpo¬ 
sition  is 


Due  to  these  phase  factors  eu!',  the  motional  superposition  is 
not  correctly  reversed,  as  shown  in  Fig.  2,  panel  4.  The  sec¬ 
ond  beam  splitter  consists  of  the  displacements  T>(  —  af  and 
T>(  —  at)  followed  by  a  second  tt/2  pulse  on  the  carrier, 
yielding 


l<A>=  ^ll)(e 


,!<^|m)  +  e~'  n^\n)). 


(45) 


In  the  second  beam  splitter  a  —  7r  pulse  again  entangles  the 
spins  with  the  motion,  as 


1^)=  (e--0|i)|O)  +  e-'^|T)|n')).  (46) 

V2 


+  e'KI~sln^(|  j)  —  e'5\]))\aAe,‘t>—  1))].  (41)  A  final  7r/2  pulse  generates  the 


state 
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FIG.  3.  Pictorial  representation  of  the  action  of  a  phase  reservoir 
on  a  Fock  state  interferometer.  Here  we  consider  the  interferometer 
that  generates  the  state  | <A)  =  ||)(|0)  + 1 2))/ \Jl.  A  ir/2  pulse  on  the 
first  blue  sideband  generates  the  state  | <f)  =  (|l)|0)  + 1 T)|  1))/ ^2 
depicted  in  panel  1,  where  the  different  spin  states  are  indicated  by 
the  shading.  A  rr  pulse  tuned  to  the  first  red  sideband  drives  the 
transition  |T)|l)— >||)|2)  and  generates  the  desired  superposition. 
The  random  motional  phase  shift  4>  in  panel  3  adds  a  relative  phase 
factor,  |i/f)  =  |  j.)(|0)  +  e-2,^|2))/\/2,  where  the  phase  is  scaled  by 
the  difference  in  indices  of  the  Fock  states  (2  in  this  case).  A  second 
7r  pulse  generates  the  state  |  <A)  =  (|  l)|0}  +  e~2'^|  |)|  1 ))/ (panel 
4).  The  result  of  the  final  ir/2  pulse  depends  on  the  phase  </>.  Ran¬ 
dom  fluctuations  in  <j>  result  in  decreased  contrast  of  the  Ramsey 
fringes. 


\t)=1-[(e-irn*-e«s-"‘»)\l)\0) 

+  (e-,^+ei('5-m^))|t)| «')].  (47) 

The  projection  onto  ||)  yields 

pl((l>)  =  \(l\^)\2  =  ^{i  +cos[  S+  (n~m  )</>]}.  (48) 

Again,  the  phase  shift  introduced  between  the  beam  split¬ 
ter  pulses  results  in  a  shift  of  the  fringes.  The  random-phase 
shifts  cf>  are  assumed  to  be  Gaussian  distributed  with  stan¬ 
dard  deviation  <x,  and  the  ensemble  average  of  P  ( <f>)  is 

Pi  =  I[l+e-(»-»)2«^nCosi5].  (49) 

The  decoherence  scales  exponentially  with  the  square  of  the 
difference  in  the  Fock  state  indices  n  —  m,  so  that  “larger” 
superpositions  lose  coherence  exponentially  faster.  Compare 
the  scaling  of  the  decay  of  the  fringe  contrast  as  a  function  of 
n  —  m  in  Eq.  (49)  with  the  decay  of  the  off-diagonal  elements 
of  the  density  matrix  pmn  given  in  Eq.  (21).  The  results  are 
equivalent.  (See  Fig.  3.) 

C.  Connection  of  averaged  interferometer  approach 
and  master  equation 

It  is  reasonable  to  ask  in  what  sense  the  interferometer 
experiments  of  Sec.  II B  directly  probe  the  relevant  off- 
diagonal  matrix  elements  as  calculated  in  Sec.  II  A.  In  the 
case  of  Fock  states  interacting  with  the  phase  reservoir  (Secs. 
II  A  2  and  II  B  3),  the  connection  is  quite  simple,  in  other 
cases,  the  connection  is  not  as  straightforward. 
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1.  Spin  coherence 

In  this  section  we  argue  that  the  contrast  of  the  interfer¬ 
ence  fringes  is  a  measure  of  the  coherence  of  the  combined 
spin  and  motional  quantum  state,  and  that  the  spin  part  can 
be  negligible  or  can  be  removed.  To  address  this  issue,  we 
consider  the  interferometer  experiment  in  the  case  of  no 
added  noise.  We  take  the  example  of  the  cat  state  interfer¬ 
ometer  discussed  in  Sec.  II  B  and  write  the  cat  state  of  the 
system  after  the  first  rr/2  and  displacement  pulses  (assumed 
to  be  short)  as  a  density  matrix  [the  wave  function  for  C 
=  1  is  given  in  Eq.  25]: 

P=  2  [|«l)l IXI K«[l  + 1  «t)ITXT|<«tI 

+  C|ai)||)(TK«Tl+C*|«T)|t)U|(ai|],  (50) 

where  C=CsCm  quantifies  the  coherence  that  has  contribu¬ 
tions  from  both  the  spin  (Cs)  and  motion  (Cm).  This  separa¬ 
tion  of  the  contrast  into  a  product  of  spin  and  motion  contri¬ 
butions  is  valid  for  reservoirs  in  which  the  spin  and  motion 
are  not  coupled,  such  as  the  amplitude,  natural,  and  phase 
reservoirs.  Moreover,  separate  experiments  verify  that  the 
spin  is  not  disturbed  by  or  during  the  time  of  application  of 
our  engineered  phase,  amplitude,  and  zero-temperature  res¬ 
ervoirs,  so  that  in  these  experiments  the  fringe  contrast  is  a 
simple  probe  of  the  change  in  motional  coherence.  The  only 
reservoir  experiments  in  which  the  spin  decoherence  is  not 
negligible  are  the  natural  reservoir  experiments,  since  for  all 
of  our  applied  reservoirs  the  time  of  application  is  short.  For 
the  natural  reservoir  experiments,  we  are  forced  to  simply 
wait  a  time  t  for  the  ambient  noise  to  act  on  the  system.  Note 
that  the  separation  of  coherences  is  valid  if  the  spin  under¬ 
goes  only  phase  decay,  with  no  random  spin-flip  transitions, 
which  is  verified  in  our  experiments.  To  within  experimental 
errors,  we  can  adequately  characterize  spin  decoherence  as 
exponential,  CJ  =  exp(  —  yst). 

In  the  second  part  of  the  cat  state  interferometer,  the  sec¬ 
ond  displacement  is  applied,  undoing  the  first  displacement, 
leaving 


p=t[|o>|I>UI<o|  +  |o)|t><TI<o| 

+c|o)U)(r|(o|+c*]o)|r)UKo|].  (51) 

The  final  7r/2 -pulse  with  phase  S  transforms  the  spins  as  in 
Eqs.  (30),  and  the  detection  of  |  J, )  with  trace  over  the  motion 
leaves: 


1  —  Re(C)cos  S 

A= - j - '  (52) 

Thus  the  contrast  of  the  Ramsey  fringes  depends  directly  on 
the  coherence. 


2.  Cats,  amplitude  noise 

The  interferometer  calculations  presented  in  Sec.  II  B  pre¬ 
dict  the  outcome  of  particular  experiments.  The  master- 
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equation  solutions  in  Sec.  II  A  do  not.  In  this  section,  we 
take  the  solution  of  the  master  equation,  in  terms  of  the  time 
evolution  of  a  density  matrix  and  derive  a  prediction  for  an 
interferometer  experiment.  The  result  is  shown  to  be  equiva¬ 
lent  to  the  average-over-classical-variables  approach  in  Sec. 
II  B.  We  do  this  for  the  Schrodinger  cat  interferometer  sub¬ 
ject  to  amplitude  noise.  The  solution  to  the  master  equation 
was  given  in  Sec.  II  A  1 .  It  is  a  straightforward  matter  to  add 
the  spin  part  to  this  solution,  which  must  be  done  for  the 
interferometer  calculation.  Starting  from  Eq.  (14),  we  divide 
the  initial  density  matrix  [written  as  an  initial  P  function 
P0(/i,v*)]  into  four  parts,  corresponding  to  the  initial  state 

l<A)  =  (ll)l«|}+lT>K))/^2: 

Po(/i,v*)  =  Pio\^i'*)  +  Po\^*)  +  Po\^i'*) 

+  P$\p,v*)  (53) 

with 

P(o\fi,v*)=  ^  <5(,u- oq)^  i'-a1)|J.)(l|,  (54) 

P{Q2\fj,,v*)=  -5(/r-aT)5(v-aT)|t)(T|,  (55) 

fo3,(  V’ v* )  =  y  ( «T  I  «|)  <5(  P-  ~  ai )  <*(  v~  «t )  1 1 )(  T I > 

(56) 

ptf\n,v*)  =  [ptf\v,v*)y-  (57) 


that  transforms  the  spins  as  in  Eq.  (30)  and  measure  the 
probability  that  the  ion  is  in  the  spin-down  state,  as  discussed 
in  the  context  of  Eq.  (23).  Thus  we  are  interested  only  in  the 
component  of  p  given  by  Pj  =  (J.]p/U),  where  pf  is  the  den¬ 
sity  matrix  after  the  final  displacement  and  Ramsey  pulses. 
Furthermore,  since  we  do  not  measure  the  motional  state,  a 
trace  over  the  motional  degrees  of  freedom  is  required: 

Bi(f)  =  Tr„pI(f)  =  2  (n\pi(t)\n).  (64) 

n 

The  quantity  p  j  is  given  by 

Pl(t)  =  p{1\t)  +  p(i\t)  +  pf\t)  +  p(l4\t )  (65) 

with 

p\l\t)=l-V(ale-ie)V(al)pT(t)V\al)VHale-ie), 

(66) 

pf\t)  =  ^p{a]e-,e)V(a])pT(t)V'<  (a^V1  (a^e-'8), 

(67) 

pf)(t)=-l-eiS{a]\al)l-[e~y'n+2N^D{a[e-ie) 

XV(al)pT(t)V\a^)V\a^-ie),  (68) 

p\4\t)  =  [pf\t)V.  (69) 


From  Eq.  (14)  this  leads  to  four  contributions  to  p(f) 

P(0  =  P(1)(0  +  P(2)W  +  P°\t)  +  P(4)W  (58) 

with 


p^\t)=l-V(al)pT(t)V\ai)\l)(i\, 

(59) 

pV\t)=1-V(a1)pT(t)VH*1)  ItXtl, 

(60) 

(aT|a1)1^[e'y'/1+2A,(r)]X>(51)pr(f)£,t( 

«T)U)(tl 

(61) 

p(4)(0  =  [p(3)(f)]t, 

(62) 

where  aj  t  are  defined  as  in  Eqs.  (15)  and  (16).  To  complete 
the  experiment,  we  apply  a  second  displacement,  which  is 
given  by  the  operator 


In  what  follows,  we  assume  that  ffj  ,  are  real.  Performing 
the  trace,  we  find  that  Tr„p(l^  =  Tr„pf2^  =  1/4  and  that 
Tr„p(3'4'  gives  rise  to  the  interference.  Using  the  identity 

V{K)V{\)  =  e^KX*-K*x)V{K+\),  (70) 


substituting  for  pT  from  Eq.  (17)  and  using  the  closure  rela¬ 
tion  Sm|m)(m|  =  1, 


T r„pf)(f) 


1 

1  +N(t)  6 


i(a^  —  aj_)(aj_  +  <x-|-)sin  6 


j  N(t)  \ 

\1  +N{t)J 


A„( \P\2)e 


l/S|2, 


(71) 


where  /3=a^e  ,e+a^  —  a<e  ,e—  a.  and  £  is  a  Laguerre 
polynomial  which  comes  from  (m\D(  K)\m) 
=  exp( — |  /c|2/2)  £m  ( |  /c  | 2) .  Using  ZmemCm(x)  =  (  1 

—  e)  'exp[  — e<7(l  —  e)]  and  further  simplifying. 


Pl(t) 


^ _ —  (q:|  —  ap2[l  +  n(l  —  e  yt)  +  e  yt/2cos  0] 

2  2 e 


l)(l\V{a[e-ie)  +  \])(]\V{^e->e),  (63) 


Xcos[<5+e  yf/2(a2  —  a2)sin  ff\.  (72) 


where  6  is  the  phase  of  the  second  displacement  with  respect  For  all  of  our  experiments  we  can  consider  the  high- 
to  that  of  the  first.  We  then  apply  the  final  tt/2  Ramsey  pulse  temperature,  short-time  limits:  n>  1  and  yt<  1 .  In  this  limit. 
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mm 


FIG.  4.  Schematic  diagram  of  the  electrodes  of  the  linear  trap 
(not  to  scale).  The  trap  electrodes  are  formed  by  evaporating  gold 
onto  alumina  substrates.  The  outer  segmented  electrodes  are  the 
endcaps,  while  the  long  unbroken  electrodes  carry  rf.  The  axial  (z) 
direction  is  parallel  to  the  rf  electrode.  The  two  separate  trap  wafers 
are  separated  by  200  /rm  with  alumina  spacers  (not  shown). 


2p2Pi/2 


L-i= IT>  1  2 

f  2s  2Si/2 

2,-1  J 

FIG.  5.  Simplified  level  diagram  of  9Be+  (not  to  scale).  The 
transitions  labeled  D  are:  Dl:  optical  pumping  (also  called  red  Dop¬ 
pler),  D2:  detection/cooling,  D3:  repumping.  The  transitions  labeled 
R  are  the  two  arms  of  the  two-photon  stimulated-Raman  transition. 
A  is  the  detuning  between  the  2 p  2P1/2  and  the  Raman  virtual  level 
and  w0/2tt=  1.25  GHz  is  the  2s  2S1/2  hyperfine  splitting.  The  fine- 
structure  splitting  of  the  P  state  is  197  GHz.  The  transition  wave¬ 
length  from  the  S  to  P  state  is  313  nm. 


Pl(t) 


(aj  — +«y?  +  cos  0] 


Xcos[<5+  (1  —  yt/2)(a2  —  a2) sin  0]}.  (73) 


Normally,  the  second  displacement  is  the  reverse  of  the  first 
{6=  if)  and  the  Ramsey  fringes  in  Eq.  (33)  are  recovered 
with  the  association  uyf— >2cr2.  The  scaling  with  the  square 
of  the  size  of  the  initial  state  is  evident.  The  constant  term  in 
the  exponential  measures  a  residual  overlap  between  the  ini¬ 
tial  state  and  the  final  state,  present  even  in  the  infinite  time 
limit. 

Similar  procedures  to  this  for  other  types  of  initial  states 
and  reservoirs  yield  results  that  agree  with  those  from  the 
interferometer  averaging  over  random  displacements/phases 
approach  given  in  Sec.  II  B. 


ITT.  EXPERIMENT 
A.  Ion  trap 

We  use  a  miniaturized  version  of  the  linear  Paul  trap  [21]. 
A  diagram  of  the  trap  structure  is  shown  in  Fig.  4.  It  consists 
of  two  alumina  wafers  separated  by  200  /zm  with  alumina 
spacers  (not  shown).  Each  wafer  has  a  slot  2  mm  long  by 
200  /iin  wide.  Gold  is  deposited  onto  the  edges  of  the  slot  to 
form  electrodes.  One  side  of  the  slot  is  the  rf  electrode,  and 
the  other  side  is  divided  into  three  segments  to  which  differ¬ 
ent  static  potentials  are  applied  to  make  a  static  well  along 
the  z  axis.  The  segmented  electrodes  are  formed  by  cutting 
two  side  slots  of  width  20  /zm.  The  side  slots  are  separated 
by  200  /zm,  forming  two  “endcaps”  and  one  “middle” 
electrode.  When  the  wafers  are  sandwiched  together  with  the 
spacers,  the  four  interior  edges  of  the  main  slots  approximate 
the  four  wires  of  a  linear  quadrupole  trap.  The  microfabri- 
cated  structure  described  here  allows  miniaturization  of  the 
linear  trap  and  correspondingly  high-trap  strength. 

The  linear  trap  is  mounted  at  the  end  of  an  rf  quarter- 
wave  coaxial  resonator  that  supplies  the  required  rf  voltage 
[22].  The  trap  and  coaxial  resonator  are  contained  within  an 
UHV  enclosure,  as  in  Ref.  [22].  The  two  rf  electrodes  are 


connected  to  the  tip  of  the  coaxial  center  conductor.  Each  of 
the  static -potential  segments  are  connected  to  UHV 
feedthroughs  through  a  pair  of  200  kHz  RC  low-pass  filters 
in  series.  The  coaxial  resonator  is  30  cm  long  and  partially 
filled  with  alumina  to  lower  its  resonance  frequency  to 
Q,tI2tt=  113  MHz.  With  8  W  of  rf  power,  the  radial  secular 
frequencies  are  approximately  13  MHz.  The  axial  potential  is 
controlled  by  the  static  potentials  applied  to  the  segmented 
electrodes.  Positive  potentials  are  applied  to  the  endcaps  to 
confine  positive  ions  in  the  middle  of  the  trap.  For  the  data 
presented  in  this  paper,  we  restrict  our  attention  to  the  axial 
mode  of  motion.  With  30  V  applied  to  the  endcap  segments 
and  the  middle  segments  and  rf  electrodes  held  at  static 
ground  potential,  the  axial  frequency  is  a>J2ir=  1 1.3  MHz. 

B.  Raman  transitions  and  cooling  the  ion 

The  level  structure  of  9Be+  is  shown  in  Fig.  5.  We  select 
two  states  from  the  ground-state  hyperfine  manifold:  the  |  F 
=  2,mF=—  2)  state,  which  we  label  |J_),  and  the  \F=l,mF 
=  —  1 )  state,  which  we  label  |  f ).  The  ion  is  Doppler  cooled 
by  applying  circularly  polarized  light  detuned  ~8  MHz  be¬ 
low  the  strongly-allowed  (17277=19  MHz)  cycling  transi¬ 
tion  2s  25’1/2|2,  —  2)— >2p  2P3/2|3,—  3)  (beam  D2  of  Fig.  5). 
Due  to  imperfect  a  polarization  the  ion  optically  pumps  to 
the  1 1,  —  1 )  state  with  a  probability  —1X10  4 .  To  prevent  a 
loss  of  cooling  from  these  decay  events,  we  apply  laser  light 
resonant  with  the  transition  2s“S1/2|l,  1)  —>2p  P 1/21 2, 

—  2) .  This  transition  is  labeled  Dl  in  Fig.  5.  After  Doppler 
cooling,  the  mean  occupation  number  of  the  axial  mode  of 
motion  is  . 

We  efficiently  detect  (Fig.  6)  the  |J.)  internal  state  of  the 
ion  by  applying  Doppler  beam  D2  with  no  Dl  light  [18].  The 
fluorescence  collection  system  has  a//1.2  aperture,  and  the 
light  is  transferred  onto  the  photocathode  of  a  photomulti¬ 
plier  tube  (PMT).  The  PMT  quantum  efficiency  is  «»30%. 
After  all  losses  at  windows  and  filters,  the  overall  scattered- 
photon  detection  efficiency  is  measured  to  be  «*1.5X10~3. 
For  detection,  the  laser  is  at  the  same  detuning  as  for  Dop¬ 
pler  cooling  and  the  intensity  is  kept  well  below  saturation 
7~0.1/sat.  During  a  200-/zs-dctection  period,  we  collect  an 
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Photons  collected  in  200  ns 


FIG.  6.  Detection  efficiency.  The  graphs  show  histograms  of 
photons  collected  in  a  200  /is  interval  for  2000  repetitions  of  the 
experiment  for  the  two  initial  states  ||)  and  ||).  Note  that  if  a 
discriminator  were  placed  at  the  bin  associated  with  n  —  3  photons, 
then,  on  a  per  experiment  basis,  the  state  of  the  ion  would  be 
distinguishable  with  a  97%  quantum  efficiency. 


average  of  10  photons  if  the  ion  is  in  the  ||)  state,  and  0.1 
photons  if  it  is  in  the  1 j )  state.  The  photon  distributions  are 
Poissonian.  By  averaging  the  number  of  photons  detected  in 
successive  repetitions  of  a  given  experiment  and  normalizing 
to  the  count  rate  when  the  ion  is  prepared  in  a  control  ex¬ 
periment  in  the  ||)  state  (by  simple  optical  pumping),  we 
determine  the  probability  P  j  . 

In  previous  experiments  using  miniature  spherical  quad- 
rupole  ion  traps  (see  e.g..  Ref.  [24]),  ion  internal  state  dis¬ 
crimination  was  limited  by  sidebands  on  the  2  s  2S,/2F=1 
— >2 p  2P|/2 ,  F=2  transition  due  to  ion  micromotion  at  fre¬ 
quency  O  /  .  Thus  a  laser  tuned  to  transition  D2  would 
weakly  drive  an  ion  in  the  | f )  state.  The  miniature  linear  trap 
addresses  the  issue  of  micromotion  in  two  ways.  First,  there 
is  negligible  micromotion  associated  with  the  confinement 
along  the  axis  of  the  trap.  Second,  the  close  proximity  of  the 
static-potential  electrodes  permits  accurate  cancellation  of 
stray  static  fields,  resulting  in  negligible  sidebands  on  the 
transitions  from  micromotion.  This  permits  state  discrimina¬ 
tion  of  97%  in  a  single  experiment.  This  is  shown  in  Fig.  6. 
Even  in  the  absence  of  micromotion,  detection  efficiency  is 
limited  by  off-resonant  optical  pumping  of  the  1 1 )  state  by 
the  detection  laser  [23]. 

We  drive  coherent  stimulated  Raman  transitions  with  a 
pair  of  laser  beams  (such  as  R1  and  R2  shown  in  Fig.  5) 
detuned  A^ir^lO  to  20  GHz  from  the  2  s  2S,1/2— >2/?  2P m 
transition  near  A  =  3 1 3  nm.  The  Raman  beams  are  derived 
from  a  single  dye  laser  beam  that  is  frequency  doubled,  and 
the  difference  frequency  between  the  two  Raman  beams  is 


generated  with  several  acousto-optic  modulators  in  series. 
We  employ  three  types  of  Raman  transitions,  distinguished 
by  beam  geometry  and  difference  frequency  [16,25].  The 
beam  geometry  sets  the  projection  of  photon  momentum 
onto  the  axis  of  the  trap  (z),  characterized  by  the  wave  vec¬ 
tor  difference  <5k  =  (k[  — k2)-z.  The  coupling  of  this  photon 
momentum  to  the  motion  of  the  ion  is  parametrized  by  the 
Lamb-Dicke  parameter  r;=  Sk  z0,  where  z0  =  \h/2mtjj~=7 
nm  is  the  harmonic-oscillator  characteristic  length  for  axial 
motion. 

The  three  types  of  Raman  transitions  are  as  follows,  (i) 
Copropagating  beams  (Sk^O)  drive  a  motion-independent 
“spin-flip”  transition,  |f)|n)<-»  |  J, )  |  n ) .  The  frequency  differ¬ 
ence  between  the  Raman  beams  is  set  equal  to  the  frequency 
difference  between  the  states,  <u0=  &>/,/+  wz,  where 
(x>hfllTT^\.25  GHz  is  the  hyperfine  splitting  and  oj//2tt 
=  12  MHz  is  the  Zeeman  splitting  due  to  an  applied  mag¬ 
netic  field  |  B|  ^  0.57  mT.  (ii)  Beams  oriented  45°  to  the  trap 
axis  and  90°  to  each  other  result  in  a  wave  vector  difference 
pointing  along  the  trap  axis  with  8k  =  2yjl  tt/X  and  couple  to 
the  axial  motion.  For  the  conditions  of  this  experiment,  77 
=  0.20.  These  beams  at  difference  frequencies  a>Q+a>An, 
couple  the  states  ||)|«)  and  |j")|n  +  An).  In  the  Lamb-Dicke 
limit,  (?72(n)<gl),  the  strength  of  the  sideband  coupling  is 
proportional  to  r)^nV  The  Raman  beams  in  cases  (i)  and  (ii) 
are  linearly  polarized,  with  one  beam  polarized  parallel  to 
the  quantization  axis  (Am  =  0)  and  the  other  polarized  per¬ 
pendicular  (A m  =  ±1).  Note  that  only  the  tr  component  of 
the  latter  beam  is  required  to  drive  the  Raman  transition,  but 
the  use  of  linear  polarized  light  reduces  the  Stark  shifts  rela¬ 
tive  to  the  case  of  using  circularly  polarized  light,  (iii)  Beams 
oriented  as  in  (ii),  with  difference  frequency  set  equal  to  the 
axial  secular  motion  frequency  oj  approximate  the  harmonic- 
oscillator  displacement  operator  V(a),  defined  by  the  rela¬ 
tion  T>{  a) 1 0)  =  |  a),  and  | «)  is  a  coherent  state.  In  the  Lamb- 
Dicke  limit,  the  displacement  |a|=  r/Cl-pt  is  proportional  to 
the  duration  of  the  laser  pulse  t  and  the  Rabi  frequency  f  i-D , 
and  0  is  set  by  the  phase  and  polarization  of  the  applied  laser 
fields.  In  this  case,  the  two  Raman  beams  are  linearly  polar¬ 
ized  such  that  the  polarization  vectors  are  mutually  orthogo¬ 
nal  to  both  each  other  and  to  the  quantization  axis.  This 
cancels  Stark  shifts  to  a  high  degree  and  results  in  the  two 
spin  states  being  simultaneously  displaced  in  opposite  direc¬ 
tions  such  that  a^  =  —  a^/ 2. 

The  ion  is  first  prepared  in  the  ground  state  of  motion 
with  resolved  sideband  Raman  cooling  [18].  Raman  cooling 
consists  of  repeatedly  driving  the  first  red-Raman  sideband 
of  the  axial-motion-sensitive  transition,  case  (ii)  with  An 
=  — 1,  followed  by  optical  pumping  to  the  |J,)  state  after 
each  Raman  pulse.  The  optical  pumping  is  achieved  by  ap¬ 
plying  beams  D1  and  D3  (see  Lig.  5).  The  rate  at  which  the 
ion  undergoes  the  red  sideband  Raman  transition  depends  on 
the  initial  state  of  motion  as  flrsb°c  \Jn.  The  duration  of  the 
successive  Raman  pulses  is  progressively  lengthened,  so  as 
to  successively  drive  77  pulses  on  the  n,-— >n,-—  l,n,—  1— >n,- 
—  2,  etc.,  transitions,  where  nt  is  the  total  number  of  cooling 
pulses.  The  final  pulse  of  this  sequence  is  selected  to  drive  a 
77  pulse  on  the  |J.)|l)— >||)|0)  transition.  After  each  Raman 
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FIG.  7.  Calibration  of  displacement  amplitude 
for  (a)  uniform  electric-field  drive  and  (b)  laser 
drive.  The  y  axis  is  in  terms  of  the  coherent-state 
amplitude  a,  determined  as  described  in  the  text. 


pulse,  we  drive  the  D 1  and  D 3  transitions  for  5  /xs,  optically 
pumping  any  population  in  the  |f)  state,  or  the  (2,-1)  state 
to  the  ||)  state.  Typically,  after  6  cycles  of  Raman  cooling, 
the  ion  has  a  probability  of  95%  to  be  in  the  ground  state 

ll>|o>. 

C.  Ion  interferometry:  techniques 

In  the  experiments  reported  here,  we  measure  the  coher¬ 
ence  of  quantum  superpositions  with  single-atom  interferom¬ 
etry,  analogous  to  that  used  in  our  previous  work  [10]  and 
discussed  in  Sec.  II  B.  The  motional  state  of  the  ion  is  split 
into  a  superposition  of  two  parts.  The  “beam  splitters”  of 
our  interferometers  consist  of  a  pair  of  stimulated  Raman 
pulses.  We  employ  two  types  of  interferometers,  both  dis¬ 
cussed  in  Sec.  II  B. 

1.  Ion  interferometry:  Cat  states 

The  Schrodinger-cat  interferometer  starts  with  the  usual 
ground-state  (||)|0))  preparation.  The  first  pulse  is  then  a 
7t/2  pulse  on  the  motion-independent  transition  (i),  with 
drive  time  typically  0.3  /xs,  which  generates  an  equal 

spin  superposition,  ||)|0)— KlD  + lt))]0)/V2  as  indicated 
schematically  in  Fig.  1,  panel  (1). 

We  drive  Raman  transition  (iii)  to  excite  the  motion  asso¬ 
ciated  with  each  spin  state  into  a  coherent  state,  resulting  in 
the  state  of  Eq.  (25).  We  call  this  the  displacement  pulse 
[Fig.  1,  panel  (2)].  We  vary  the  length  of  the  Raman  drive 
pulse  in  order  to  vary  |  a ,  |  and  « ■  | .  We  independently  mea¬ 
sure  |  or  1 1  and  |aq|  in  the  following  manner  [17]  [which  is 
also  discussed  in  the  context  of  Eqs.  (85)  and  (86),  below]. 
The  ion  is  prepared  in  a  coherent  state,  and  then  Raman 
transition  (ii),  tuned  to  the  first  blue  sideband,  is  driven  for  a 
variable  time  tp.  After  this  “probe”  Raman  transition,  we 
determine  the  probability  P \ftp).  The  resulting  curve  may  be 
decomposed  into  sums  of  sinusoids  that  correspond  to  An 
=  1  transitions  between  Fock  states.  The  Fock  state  distribu¬ 


tion  is  then  fit  to  that  of  a  coherent  state,  thereby  extracting 
|a|  [17].  The  values  of  \a  \  and  |at|  are  measured  sepa¬ 
rately,  by  preparing  the  ion  in  either  the  1 1)|0)  state  or  in  the 
IT ) 1 0)  state,  applying  the  coherent  drive  pulse,  and  then 
probing  and  detecting  as  just  described.  Curves  of  |  aq  J  vs 
the  length  t  of  the  displacement  pulse  are  shown  in  Fig.  7(b). 
The  ratio  |aq  |/|aq|  =  2.0±0.1,  averaged  over  the  data  set,  is 
consistent  with  the  expected  value  of  2.  The  slope 
d\Aa\/dt=  r] LI  is  consistent  with  the  expected  value. 

The  Schrodinger-cat  interferometer  is  completed  by  ap¬ 
plying  a  second  coherent  displacement  pulse  180°  out  of 
phase  from  the  first  displacement  pulse  [Fig.  1,  panel  (4)]. 
The  phase  of  the  displacement  is  set  by  the  phase  of  an 
oscillator  that  drives  an  acousto-optic  modulator  that  is  in 
one  of  the  Raman  beams.  All  phases  are  referenced  to  a 
master  oscillator  by  phase  locking  to  that  oscillator.  If  the 
motion  was  unperturbed  during  the  time  between  the  dis¬ 
placement  pulses,  then  the  inverse  displacement  would  re¬ 
store  the  motion  associated  with  each  spin  state  to  the  ground 
state  of  motion.  The  wave  function  would  return  to  |  tp) 
=  (U>  +  IT»|0>/V2.  A  second  tt/2  pulse  on  the  motion- 
independent  spin-flip  transition,  with  phase  S  relative  to  the 
first  7t/2  pulse,  interferes  with  the  spin  components.  The 
probability  P  ( 8)  varies  sinusoidally  with  S.  To  record 
P l(S),  we  step  the  frequency  of  the  oscillator  controlling  the 
7t/2  pulses  and  bin  the  fluorescence  counts  at  each  step. 
Here,  S  is  the  accumulated  phase  difference  caused  by  the 
frequency  difference  between  the  oscillator  frequency  v  and 
v{]  as  in  normal  Ramsey  spectroscopy.  We  typically  repeat 
the  experiment  100  times  at  each  step.  Note  that  the  resulting 
fringe  spacing  is  much  smaller  than  the  width  of  the  Rabi 
pedestal,  meaning  that  the  predominant  effect  is  a  change  in 
phase,  without  introducing  significant  inaccuracy  in  the  7t/2 
pulse.  We  typically  record  three  oscillations  of  P ,  .  The 
fringes  are  fit  with  a  sinusoid  A+B  cos(&>—  co0)T,  where  T  is 
the  time  between  spin  pulses.  The  contrast  of  the  fringes  is 
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defined  as  BIA .  For  the  basic  interferometer  just  described, 
with  no  perturbations  of  the  motion  between  the  “beam 
splitter”  pulses,  the  contrast  is  typically  0.6-0. 8.  This  con¬ 
trast  differs  from  unity  for  several  reasons.  (1)  Errors  in  the 
7t/2  spin-flip  pulses,  primarily  due  to  laser  power  fluctua¬ 
tions.  (2)  Fluctuating  magnetic  fields  that  reduce  the  spin 
coherence  through  fluctuating  Zeeman  shifts.  (3)  Motional 
heating  during  and  between  the  displacement  pulses,  with 
further  small  contributions  due  to  inaccuracies  in  reversing 
the  initial  displacement  pulse.  These  effects  are  expected  to 
be  random  and  should  not  influence  the  additional  decoher¬ 
ence  caused  by  the  applied  reservoir  (see  Sec.  II  C). 

2.  Interferometry:  Fock  state  interferometer 

The  Fock  state  interferometer  beam  splitters  combine  a 
7t/2  pulse  with  a  7r  pulse  on  transitions  of  either  type  (i)  or 
(ii),  as  shown  in  Table  I.  As  an  example,  consider  the  gen¬ 
eration  of  the  superposition  1 1)(|0)  +  \2))/\j2.  A  Tr/2-pulse 
on  the  first  blue  sideband,  transition  type  (ii)  with  An  =  +  l 
and  duration  7’7I./2S=  1  /x s,  drives  the  transition  ||)|0) 
-(U>|0>  +  |T>|1»/ a/2.  A  7 t  pulse  on  the  first  red  sideband, 
transition  type  (ii)  with  A«  =  — 1  and  duration  '/’_«=  1 .5  /xs, 
drives  the  transition  |T)|  1)— >||)|2).  This  second  pulse  does 
not  affect  the  population  in  the  ||)|0)  state.  The  result  of 
these  two  pulses  is  the  desired  state  |  *A)=  1 |)(|0)  + 1 2 ) )/  \/2. 
The  second  beam  splitter  recombines  the  two  parts  of  the 
wave  function  by  driving  a  second  t t  pulse  on  the  red  side¬ 
band  followed  by  a  n/2  pulse  on  the  blue  sideband.  By 
sweeping  the  frequency  of  the  blue  sideband  oscillator  we 
generate  Ramsey  fringes  in  the  probability  P  ,  .  Contrast 
(without  coupling  to  the  reservoir)  is  limited  by  the  same 
mechanisms  as  in  the  Schrodinger-cat  interferometer. 

D.  Engineering  reservoirs:  techniques 
1.  High-temperature  amplitude  reservoir 

As  discussed  in  the  Introduction,  the  motion  of  a  trapped 
ion  couples  to  environmental  electric  fields  E  through  the 
potential  U=  —  c/x-  E,  where  x  is  the  position  of  the  ion  rela¬ 
tive  to  its  equilibrium  position  (proportional  to  the  amplitude 
of  motion)  and  q  is  the  charge  of  the  ion.  If  E  is  due  to  a 
reservoir  of  fluctuating  field  modes,  then  E°c2, €,(£,  +  £>/), 
where  h , ,  b\  are  the  lowering  and  raising  operators  of  the 
field  modes  (as  in  Sec.  II  A  1).  A  classical  coherent  drive  (for 
which  E  is  a  narrow-band  sine-wave)  applied  to  the  trap 
electrodes  at  the  ion  axial  motional  frequency  a>  results  in  a 
displacement  of  the  motional  state  proportional  to  the  size 
and  duration  of  the  applied  field.  This  is  shown  in  Fig.  7(a) 
for  two  different  durations  of  the  applied  field.  We  simulate 
the  effects  of  the  high-temperature  amplitude  reservoir  dis¬ 
cussed  in  Sec.  II A  1  by  applying  random  electric  fields 
along  the  axis  of  the  trap,  whose  spectrum  is  centered  on  the 
ion  axial  motion  frequency  a>. 

To  generate  the  required  axial  fields,  we  apply  a  potential 
to  one  of  the  endcap  electrodes.  This  generates  both  an  axial 
field  and  a  small  radial  field.  We  ignore  the  effect  of  the 
radial  field  since  we  are  insensitive  to  motion  in  the  radial 
direction.  This  field  also  modulates  the  axial  potential  and 


gives  rise  to  an  accompanying  (random)  phase  shift;  how¬ 
ever,  for  the  experimental  conditions  this  phase  shift  is  neg¬ 
ligible.  We  start  with  a  10  MHz  bandwidth  noise  source  that 
rolls  off  at  6  dB/octave  starting  at  10  MHz.  We  filter  the 
noise  such  that  the  resulting  noise  spectrum  has  3  dB  points 
at  9  MHz  and  11.5  MHz,  and  the  noise  rolls  off  at  6  dB/ 
octave  above  and  below  these  points.  This  spectral  noise  is 
applied  to  the  trap  electrodes  through  a  network  of  RC  low- 
pass  filters  designed  to  prevent  such  fields  from  reaching  the 
static  trap  electrodes.  The  filters  effectively  limit  ambient 
noise  from  reaching  the  electrodes,  but  large  deliberately  ap¬ 
plied  fields  can  still  be  effective.  The  noise  is  applied  for 
3  /xs  between  the  Ramsey  zones.  The  total  time  for  the  Ram¬ 
sey  experiment  (the  time  between  the  tt/2  spin-pulses)  is 
about  T= 20  /xs  and  is  held  fixed.  For  these  filters,  trap  ge¬ 
ometry  and  application  time,  it  can  be  shown  that  this  type  of 
drive  correctly  simulates  an  infinite-bandwidth,  amplitude- 
only  reservoir. 

2.  Phase  reservoir 

A  phase  reservoir,  as  described  in  Sec.  II  A  2,  is  simulated 
by  random  variations  in  the  trap  frequency  oj,  changing  the 
phase  of  the  ion  oscillation  without  changing  its  energy.  We 
realize  this  coupling  experimentally  by  adiabatically  modu¬ 
lating  the  trap  frequency.  The  random  potential  fluctuations 
are  filtered  through  a  low-pass  network  with  a  cutoff  fre¬ 
quency  well  below  co  to  maintain  adiabaticity  and  to  avoid 
any  power  at  the  trap  secular  frequency,  which  may  cause 
amplitude  noise.  The  potential  fluctuations  are  applied  to  one 
of  the  middle  electrodes,  which  is  symmetric  with  respect  to 
the  axial  direction,  so  as  to  produce  field  gradients  and  mini¬ 
mize  any  applied  uniform  axial  fields  further  reducing  the 
possibility  of  changing  the  motional  energy.  Applying  a  po¬ 
tential  to  just  one  middle  electrode  does  result  in  a  radial 
field,  but  as  before,  we  isolate  the  axial  motion  and  ignore 
any  radial  motion. 

In  the  phase  reservoir,  the  ion  follows  the  change  in  trap 
frequency  in  the  sense  that  there  is  no  change  in  the  energy 
of  the  trap  motion,  however  the  phase  is  either  advanced  or 
retarded  relative  to  the  unperturbed  motion.  We  start  with  a 
10  MHz  pseudorandom  noise  source  and  filter  the  noise,  cre¬ 
ating  a  spectrum  that  is  flat  between  1  kHz  and  100  kHz  with 
a  roll-off  of  12  dB/octave  above  100  kHz  and  6  dB/octave 
below  1  kHz.  The  phase  deviation  of  the  ion  motion  is  re¬ 
lated  to  the  voltage  deviation  on  the  electrodes.  The  trap 
axial  frequency  is  related  to  the  voltage  difference  between 
the  endcap  and  middle  segments  by  the  relation  a>,oc^y 0. 
The  deviation  of  the  trap  frequency  is  given  by  Sa>/coz 
=»  (5V74l/o  since  only  one  electrode  is  driven.  The  phase  shift 
of  the  ion  motion  is  the  integral  of  the  frequency  deviation 
over  the  time  of  the  noise  pulse  Tn ,  Stf>=  J "  8a>:dt 
=  fTQna>z(8VIAVQ)dt.  Typically  Tn=  10  /xs.  We  cannot  di¬ 
rectly  measure  the  voltage  on  the  trap  electrodes  since  some 
of  the  filters  between  our  noise  generator  and  the  electrodes 
are  within  the  vacuum  enclosure.  However,  on  a  mock  up  of 
the  trap  electrodes  and  the  associated  filters,  we  measured  the 
voltage  deviations  with  an  integrating  circuit.  These  mea¬ 
surements  verified  that  the  phase  deviations  are  Gaussian 
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FIG.  8.  (Color)  Heating  of  the 
motional  state.  At  f  =  0  the  ion  is 
prepared  in  a  state  near  the  ground 
state  |0).  As  time  passes,  the  ion 
heats  out  of  the  ground  state,  as 
measured  by  the  energy  in  the  mo¬ 
tional  state,  which  is  proportional 
to  the  average  occupation  number 
<«>• 


t  [ms] 


distributed,  with  a  width  that  varies  with  the  applied  noise 
voltage,  and  that  negligible  voltage  noise  is  added  at  the  trap 
secular  frequency.  In  separate  experiments,  we  verified  that 
the  change  in  motional  energy  was  negligible  during  the  ap¬ 
plication  of  the  phase  reservoir. 

3.  Natural  reservoir 

The  third  type  of  reservoir  is  not  engineered,  but  is  the 
ambient,  or  “natural”  noise  in  our  trap.  We  have  previously 
observed  heating  in  miniaturized  Paul  traps  [26].  We  have 
further  determined  that  this  ambient  noise  is  dominated  by 
uniform  field  fluctuations  rather  than  field-gradient  fluctua¬ 
tions  [27].  We  do  not  currently  have  an  explanation  of  the 
heating  mechanism.  However,  thermal  or  Johnson  noise, 
background  gas  collisions,  and  several  other  sources  are  too 
small  to  explain  the  heating  rate.  We  do  know  that  this  am¬ 
bient  noise  source  adds  energy  to  the  ion  (amplitude  noise), 
and  that  the  state  of  the  ion  evolves  into  a  thermal  state  with 
steadily  increasing  temperature.  The  suspected  cause  is  fluc¬ 
tuating  patch  potentials  at  the  trap  electrodes,  a  problem  that 
is  exacerbated  by  the  small  size  of  our  trap  [26].  To  observe 
the  effects  of  this  heating  on  quantum  coherences,  we  simply 
run  the  Ramsey  interferometer  experiment  with  varying  time 
between  the  Ramsey  pulses. 

We  can  also  observe  the  effect  of  the  heating  on  the  en¬ 
ergy  of  the  motional  state.  We  find  that  due  to  our  unknown 
environmental  source  of  noisy  electric  fields,  the  ion  heats 
from  the  ground  state  at  a  rate  linear  in  time  with  d(n)/dt 
■ —  5.7 ±  1  quanta  per  ms.  Shown  in  Fig.  8  are  data  for  the 
mean  quantum  number  in  the  motion  as  a  function  of  delay 
time.  That  is,  the  ground  state  is  prepared  and  the  average 
quantum  number  is  measured  after  a  fixed  delay  time  t.  At 
each  point  the  motional  state  is  thermal  as  determined  by  the 
method  of  Refs.  [17]  and  [26]. 

4.  Zero-temperature  reservoir 

The  fourth  type  of  engineered  reservoir  was  a  zero- 
temperature  reservoir  realized  by  laser  cooling.  Some  theo¬ 
retical  results  for  this  reservoir  were  mentioned  briefly  in 


Sec.  II  A  la.  The  experimental  scheme  is  outlined  in  Fig.  9 
and  follows  a  proposal  by  Poyatos  et  al.  [28].  The  goal  is  to 
simulate  decay  of  the  states  ||)|«)  into  a  T=0  reservoir.  To 
do  this,  we  simultaneously  drive  on  the  first  red  sideband 
transitions  ||)|n)— >|"[)|«  — 1)  at  a  Rabi  frequency  Llrsb  and 
the  red-Doppler  (D1  from  Fig.  5)  from  |  f)  to  the  Pl/2  level  at 
a  rate  LlD(<sr).  This  is  effectively  a  continuous  Raman 


FIG.  9.  Engineering  a  1=0  reservoir  with  laser  cooling:  sche¬ 
matic  diagram  of  the  levels.  The  top  panel  shows  the  full  three- 
(intemal)  level  system  (dressed  by  the  motional  states),  while  the 
bottom  panel  shows  the  equivalent  two-level  system.  The  red  Dop¬ 
pler  beam  with  resonant  Rabi  frequency  flD  is  used  to  give  the  |  [ ) 
state  an  effective  linewidth  y=  .  Simultaneously  the  red 

sideband  (rsb)  is  driven  coherently  at  rate  Hrsb  to  drive  population 
out  of  the  ||)  state.  The  sum  of  the  two  applied  fields  simulates  a 
variable  bandwidth  zero  temperature  reservoir. 
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cooling  process  (in  contrast  to  the  stepwise  Raman  cooling 
process  described  at  the  end  of  Sec.  Ill  B)  which  eventually 
relaxes  the  system  to  a  state  in  which  all  population  is  in  the 
ground  state  |0)||)  [29,30].  The  |f)  level  is  given  an  effec¬ 
tive  decay  rate  ye(T~  Li  2}l F .  We  begin  by  preparing  an  initial 
state  ||)(|0)  +  |2))/\/2.  The  point  of  the  experiment  is  to 
follow  the  time  evolution  of  the  coherence  of  this  state 
(p0i  2i  as  we  label  the  relevant  off-diagonal  element  of  the 
density  matrix)  as  it  decoheres  during  the  cooling  process. 
To  gain  a  qualitative  understanding  of  the  experiment  we  can 
consider  a  reduced  problem,  as  follows. 

We  are  simulating  a  T=  0  reservoir  by  simultaneously 
driving  (resonantly)  on  the  red  sideband  and  red-Doppler 
transitions.  The  red-Doppler  drive  gives  the  |f)  state  an  ef¬ 
fective  linewidth  yeff  (see  Fig.  9).  The  system  density  matrix 
obeys  the  following  time  evolution: 

Teff 

p=i\p,H]  +  ~y [2o-_po-+-o-+o-_p-po-+o-_]  (74) 


Consider  two  limiting  cases.  In  the  first  g^Yeff,  so  that  \± 
have  an  imaginary  part.  Here  the  time  evolution  of  the  co¬ 
herence  is  an  exponentially  decaying  cosine,  reflecting  the 
fact  that  the  coherent  drive  is  faster  than  the  decay  rate: 

e  -  Teff'A  _ 

PmLOiW-*- — 2 - cos(ylmgt).  (82) 

In  the  second  case  g^ye ff,  so  that  \±  are  strictly  real.  In 
this  case, 

1  I  2s2  \ 

Pmi,0i0)^~exp  \-m - f  ,  (83) 

2  \  Yeff  / 

and  we  have  simple  exponential  decay  of  the  off-diagonal 
density-matrix  element.  It  is  instructive  to  compare  this  to 
the  T=  0  limit  of  the  amplitude  reservoir  presented  in  Sec. 
II  A  la.  From  Eq.  (11),  with  pom(0)=l/2 


with  interaction-picture  Hamiltonian 

H  =  g(a^  cr_+acr+),  (75) 


1  /  Y  \ 

Pom  =  2exP|  ~m2t  r 


(84) 


where  g  =  nrsb  is  the  Rabi  rate  of  the  red  sideband  drive 
(assumed  real)  and  cr+  (<r_)  is  the  raising  (lowering)  opera¬ 
tor  of  the  atom.  Here  we  have  neglected  the  effects  of  recoil 
heating,  valid  for  This  is  simply  the  master  equation 

for  a  driven,  decaying  two-state  atom  dressed  by  coupling  to 
a  harmonic-oscillator  system  [familiar  to  cavity  QED  as  the 
Jaynes-Cummings  model  [31]  with  atomic  (but  no  cavity) 
decay].  We  consider  the  experiment  in  which  we  make  an 
initial  state 


‘Ao  =  ll>(|0>  +  |>»»/V2  (76) 

and  we  measure  the  coherence  that  we  call  pm\m  ■  Equation 
(74)  leads  to: 

Pmi,oi=  ~i  VwgPm-iT.Oi  >  (77) 


Pm  —IT  ,0]. 


.  / —  Yeff 

i'!mgPm  i.oi  —  -7T 


■bn  — IT, 01  • 


(78) 


The  solution  to  this  set  of  differential  equations  [with  initial 
conditions  specified  by  Eq.  (76)]  is 


Pm- it,  or 


i  \jmg 


2(X+  —  X_) 


[e 


A-4 -t 


]. 


(79) 


With  the  identification  y/2^2g2/ yeff ,  Eqs.  (83)  and  (84)  are 
identical.  This  is  the  anticipated  result,  since  the  effective 
linewidth  of  the  1 1 )  state  is  simply  the  square  of  the  Rabi 
frequency  divided  by  the  bare  decay  rate,  and  in  the  case  of 
the  red  sideband  the  Rabi  frequency  scales  as  g  \in. 


IV.  DATA 

To  reiterate,  the  experiments  consist  of  a  Ramsey  type 
interferometer  experiment  in  which  the  ion  is  coupled  to  one 
of  the  engineered  reservoirs  between  the  initial  and  final 
“beam  splitter”  pulses.  The  experiment  proceeds  as  follows. 
The  ion  is  first  Doppler  cooled  and  then  Raman  cooled,  pre¬ 
paring  the  ion  in  the  ||)|0)  state.  The  state  preparation 
pulses  are  then  applied  to  generate  the  desired  superposition 
(first  beam  splitter).  A  period  of  time  elapses  during  which 
the  ion  is  coupled  to  the  reservoir.  After  this  time,  we  reverse 
the  sequence  of  pulses  used  in  state  preparation  (second 
beam  splitter).  The  resulting  internal  state  of  the  ion  depends 
on  the  frequency  and  phase  of  the  oscillators  used  in  the 
preparation  and  reversal  pulse  sequences.  The  ion  internal 
state  is  detected  as  described  in  Sec.  Ill  B.  The  contrast  of 
the  resulting  Ramsey  fringes  is  a  measure  of  the  motional 
state  coherence,  which  decreases  with  increasing  interaction 
with  the  reservoir,  either  in  time  or  in  strength  (applied  volt¬ 
age  noise). 


1 


Pm  i,0J.= 


with 


2(X  +  — X-) 


(v  ,  Yeff  \  x  f 

[v  +H>-( 

) 

ir+  2  )e 

l  2  r 

(80) 


x+  = - —  ± 


Yeff  .  1  /Yeff 


4  ~2  v  4 


-4  mg2.  (81) 


A.  Schrodinger-cat  decoherence 
1.  Cat  states,  applied  amplitude  noise 

The  effect  of  coupling  the  Schrodinger-cat  interferometer 
to  the  engineered  amplitude  reservoir  is  shown  in  Fig.  10. 
The  reservoir  is  simulated  by  driving  particular  trap  elec¬ 
trodes  with  a  noisy  potential  V  near  the  trap  secular  fre¬ 
quency  as  discussed  in  Sec.  Ill  D  1.  The  amplitude  reservoir 
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FIG.  10.  Decay  of  Schrodinger-cat  state  coherences  with  applied 
amplitude  noise.  Note  the  universal  scaling  of  the  states  to  an  ex¬ 
ponential  with  decay  constant  proportional  to  |A<?|2.  The  reservoir 
was  applied  for  3  /jls. 


is  coupled  to  the  ion  for  3  /is  between  the  beam  splitter 
pulses  of  the  interferometer.  The  variance  cr2  of  (i  (Sec. 
II  B  2a)  is  proportional  to  the  mean-squared  voltage  noise 
(V2).  Figure  10  shows  a  plot  of  the  interference  fringe  con¬ 
trast  as  a  function  of  the  applied  mean-squared  voltage, 
scaled  by  the  squared  “size”  of  the  cat  state  |Aaj2.  Decay 
curves  were  recorded  for  a  variety  of  superposition  sizes 
|Aa|,  and  all  the  data  agree  with  a  single  exponential,  as  in 
Eq.  (33).  The  initial  contrast  for  each  value  of  |Aa|  is  nor¬ 
malized  to  unity  at  ( V/2)  =  0. 

2.  Cat  states,  natural  noise 

The  natural  noise  that  gives  rise  to  the  heating  described 
in  Sec.  Ill  D  3  causes  decoherence.  For  this  experiment,  we 
simply  wait  for  some  amount  of  time  between  the  Ramsey 
zones  and  measure  the  fringe  contrast  as  a  function  of  this 
time.  The  results  are  shown  in  Fig.  11.  As  in  the  previous 
case,  the  results  are  consistent  with  a  single  exponential.  This 
time,  the  decay  constant  has  meaning  relative  to  another  ex¬ 
perimental  quantity:  it  should  be  the  same  as  that  measured 
for  the  heating.  For  example,  the  mean  quantum  number  in 
the  system  is  given  by  Eq.  (8),  from  which  y  is  interpreted  as 
the  heating  rate  of  Sec.  Ill  D  3.  The  decay  constant  relevant 
to  the  decoherence  measured  in  Fig.  1 1  is  given  by  Eq.  (73), 
with  the  same  coupling  constant  y  appearing.  The  measured 
numbers  for  y  agree  reasonably  well.  We  find  a  decay  rate  of 
yl2rr=  5.7  ±  1  quanta/ms  from  the  heating  data  (Fig.  8)  and 
yl2Tr=l .5 ± 0.7  quanta/ms  from  the  decoherence  data  (Fig. 
11).  Discrepancies  can  easily  arise  due  to  variations  in  heat¬ 
ing  rate  from  day  to  day  and  other  systematic  errors  that  are 
not  included  in  the  error  quoted  on  y,  such  as  miscalibration 
of  the  size  of  the  coherent  states  A  a.  In  addition,  there  is 
some  spin  decoherence  during  the  80  /jls  of  Ramsey  time 
used  for  the  longest  data  points  (see  the  discussion  of  Sec. 
II  C).  The  contrast  is  normalized  to  unity  at  t  =  0. 
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FIG.  1 1 .  Decay  of  Schrodinger-cat  state  coherences  in  natural 
amplitude  noise.  Note  again  the  universal  scaling  of  the  states  to  an 
exponential  with  decay  constant  proportional  to  |Aoj2,  and  that  the 
decay  constant  is  the  same  as  that  derived  from  the  heating  data  of 
Fig.  8.  This  is  discussed  further  in  the  text. 

3.  Cat  states,  phase  noise 

Results  from  coupling  of  the  Schrodinger-cat  interferom¬ 
eter  to  the  engineered  phase  reservoir  are  shown  in  Fig.  12. 
The  theory  has  only  one  free  parameter,  which  is  a  scaling  to 
convert  voltage  noise  to  phase  deviation  (see  below  for  a 
discussion  of  this).  The  theory  curves  are  from  Eq.  (44),  all 
with  the  same  scaling  parameter.  Note  that  there  is  no  simple 
universal  scaling  law  for  the  functional  form  of  the  decoher¬ 
ence  as  there  is  in  the  case  of  Schrodinger  cats  subject  to 
amplitude  noise.  The  contrasts  at  (V2)  =  0  are  normalized  to 
unity.  The  phase  noise  reservoir  is  applied  for  —10  /jls. 


FIG.  12.  Decay  of  Schrodinger-cat  state  coherences  in  an  ap¬ 
plied  phase  noise  reservoir.  Several  sizes  of  cats  are  shown.  The 
reservoir  was  applied  for  20  /is. 
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FIG.  13.  Decay  of  Fock  state  coherences  in  applied  phase  noise. 
Note  the  universal  scaling  of  the  states  to  an  exponential  with  decay 
constant  proportional  to  |  An|2.  The  reservoir  was  applied  for  10  /zs. 

B.  Decoherence  of  superpositions  of  Fock  states 
and  decay  of  Fock  states 

1 .  Fock  states,  phase  noise 

Coupling  the  Fock  state  interferometer  to  phase  noise  re¬ 
sults  in  the  simple  behavior  predicted  by  Eq.  (21).  This  is 
confirmed  in  the  data  of  Fig.  13  in  which  the  fringe  contrast 
of  the  interferometer  is  measured  for  different  values  of  the 
applied  phase  noise  reservoir.  The  data  for  different  initial 
Fock  state  superpositions  are  scaled  by  the  squared  “size” 
(An)2  of  the  superposition  in  order  to  show  the  universal 
scaling  law  of  Eq.  (21).  The  scaling  of  the  x  axis  of  the  graph 
is  a  free  parameter  in  the  theory.  In  principle,  the  size  of  the 
noise  measured  in  radians  could  be  known  from  the  variance 
of  the  applied  noise  (V2).  However,  in  practice,  unknown 
geometrical  factors  in  converting  voltage  on  a  subset  of  the 
trap  electrodes  to  change  in  trap  secular  frequency,  prevents 
a  direct  comparison.  A  rough  estimate,  nonetheless,  gives  a 
correspondence  of  about  it  radians  of  phase  noise  at  0.1  V2, 
which  is  in  rough  agreement  with  the  fit  parameters  from  the 
data.  The  voltage  noise  was  applied  for  a  10  /is  interval. 

2.  Fock  states,  amplitude  noise:  decoherence  and  decay 

The  coupling  of  superpositions  of  two  Fock  states  to  an 
amplitude  reservoir  is  particularly  relevant  to  quantum  logic 
with  trapped  ions:  Fock  state  superpositions  are  generated 
during  quantum  logic  gates,  and  the  natural  decoherence  is 
due  to  coupling  to  a  hot  thermal  amplitude  reservoir.  Deco¬ 
herence  of  superpositions  of  Fock  states  and  heating  (decay) 
of  Fock  states  limit  the  fidelity  of  quantum  operations. 

a.  Decoherence.  To  simplify  analysis,  we  restricted  our 
attention  to  superpositions  of  the  form  (\0)  +  \n))/ ^2  in 
which  Eq.  (38)  in  Sec.  II  B  2b  applies  to  a  good  approxima¬ 
tion.  Data  for  n  =  1,2,  and  3  are  shown  in  Fig.  14,  where  the 
amplitude  reservoir  is  the  natural  noise  on  the  trap  electrodes 
discussed  in  Sec.  Ill  D  3.  Using  Eq.  (38)  the  fits  Fn  for  states 
( 1 0)  + 1 « ))/ \/2  are  given  by  F„  =  exp(  —  y/)/( l+s„/)1+''  with 


FIG.  14.  Decay  of  Fock  state  coherences  in  natural  amplitude 
noise.  The  theory  curves  are  to  the  functional  form  of  Eq.  (38)  with 
an  additional  exponential  decay  to  account  for  spin  decoherence. 
The  horizontal  scaling  is  forced  to  be  the  same  for  the  three  states, 
and  the  data  are  normalized  to  unit  contrast  at  t=  0. 

sn  the  scaling  factor  forced  to  be  the  same  for  all  n.  The 
exponential  decay  accounts  for  the  loss  of  spin  coherence 
due,  for  example,  to  fluctuating  magnetic  fields,  with  ys/2ir 
~0.01  /zs  1 .  The  best  overall  fit  to  the  data  gives  ,v„  =  5/ms, 
which  is  consistent  with  the  decay  rates  quoted  in  Sec. 
IV  A  2. 

b.  Decay  of  populations;  purity  and  entropy.  In  order  to 
study  this  case  more  completely,  we  performed  measure¬ 
ments  directly  on  the  populations  (diagonal  elements  of  the 
system  density  matrix).  Starting  from  a  Fock  state  | n),  we 
applied  the  amplitude  reservoir  and  observed  the  evolution 
of  the  system.  The  procedure  was  as  follows.  We  generated 
the  state  | n)  [17],  waited  a  time  f,  during  which  the  Fock 
state  | n)  would  evolve  to  a  distribution  over  Fock  states 
{)«')}.  We  then  drove  transitions  on  the  blue  sideband  (bsb) 
for  various  Raman  probe  times  tp .  Following  this,  P  j  was 
measured  from  the  fluorescence  signal.  The  resulting  signal 
has  contributions  from  all  possible  transitions  ||)|n') 
— >|  t)|w '  +  1 ),  each  of  which  is  at  a  well-defined  frequency. 
The  bsb  curve  thus  generated  is  described  by: 

Pl(tp)=  f  +  lL  jP„[cos(mjtp)e-^b]  (85) 
with  the  Rabi  frequency  given  by 

Cll=\{n\e-i^a+a^\nF\)\  =  VLQe-,)ll2{n+\)-mVC\{ri2), 

(86) 

where  C  are  Laguerre  polynomials.  The  phenomenological 
decay  rate  yn=  —  (n+  l)0,7 To  is  used  to  model  decay  in  the 
data  due  to  Raman  beam  intensity  fluctuations  and  other 
sources  [17].  We  perform  a  singular  value  decomposition 
(svd)  [32]  on  the  data  to  extract  the  populations  P„  of  the 
various  |n)-states  at  the  end  of  the  reservoir  application  time 
t.  The  time  evolution  of  the  populations  is  shown  in  Fig.  15 
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(V2)  [volts2] 


FIG.  15.  (Color)  Decay  of  Fock  states  in  a  hot  amplitude  reservoir.  Each  graph  shows  the  time  evolution  of  the  populations  Pn  in  state 
| n)  of  the  motional  state  for  an  experiment  starting  in  the  initial  state  i/'(0)  indicated.  [The  initial  states  are  not  perfectly  generated,  as  can 
be  seen  by  P„(0).]  The  populations  are  given  in  Eq.  (10),  from  which  the  solid  lines  are  derived. 


for  different  initial  starting  states  if/(0)  =  \n0).  The  theory 
curves  are  given  in  Eq.  (10)  which  describe  the  evolution  of 
the  diagonals  of  the  density  matrix.  From  these  data,  the 
expected  qualitative  effects  are  evident.  The  theory  for  all  the 
data  of  Fig.  15  has  only  a  single  free  parameter,  which  is  the 
scaling  of  the  x  axis.  Additionally,  we  use  the  svd  popula¬ 
tions  from  the  initial  experimental  state  for  the  initial  state  of 
the  calculations  [p„„(0)  from  Eq.  (10)]. 

An  interesting  way  to  further  quantify  the  progression  to 


thermal  equilibrium  of  the  system  is  via  the  purity  and  en¬ 
tropy  of  the  (mixed)  states.  The  time  evolution  of  the  purity 
is  shown  in  Fig.  16  and  of  the  entropy  in  Fig.  17.  It  should  be 
noted  that  we  use  the  terms  entropy  and  purity  in  an  approxi¬ 
mate  sense:  the  svd  is  in  practice  only  useful  to  n^4,  so  we 
assign  a  cutoff  to  the  calculations  on  the  data.  Additionally, 
only  the  diagonal  elements  Pn=pnn  are  used.  This  is  valid 
for  initial  Fock  states  and  thermal  states  (to  which  the  ground 
state  evolves  under  the  amplitude  reservoir)  since  the  coher- 


053807-17 


Q.  A.  TURCHETTE  el  al. 


PHYSICAL  REVIEW  A  62  053807 


FIG.  16.  Evolution  of  purity  for  initial  Fock 
states  in  hot  amplitude  reservoir.  Each  graph 
shows  the  time  evolution  of  the  purity  of  the  mo¬ 
tional  state  for  an  experiment  starting  in  the  ini¬ 
tial  state  indicated.  The  purity  is  defined  in  Eq. 
(87),  from  which  the  solid  curves  are  derived. 
The  dotted  lines  are  an  untruncated  theory, 
whereas  the  solid  lines  are  truncated  at  a  maxi¬ 
mum  value  of  n  =  nmax. 


ences  are  absent.  However,  imperfections  in  state  preparation 
could  lead  to  unwanted  coherences  in  the  experimental 
states.  In  practice,  we  assume  that  such  effects  are  small,  but 
we  have  not  experimentally  verified  this  assumption.  Some 
of  the  discrepancies  between  theory  and  data  may  result  from 
this.  The  definitions  under  these  constraints  are,  for  the  pu¬ 
rity, 

n  max 

P(f)-Trp2(f)-S  P2n(t)  (87) 

n  =  0 

and  for  the  entropy, 

wmax 

5(f)-Trp(f)lnp(f)-E  P„(t)\nPn(t),  (88) 

n  =  0 

where  n  max  is  the  value  of  n  at  which  the  series  is  truncated. 
As  stated  above,  since  preparation  of  individual  Fock  states 
is  not  perfect,  it  is  necessary  to  take  into  account  the  imper¬ 
fections  in  the  initial-state  populations,  especially  as  the  ini¬ 
tial  state  gets  large.  To  match  the  data  to  the  theory,  the 
initial  diagonal  elements  of  the  density  matrix  are  measured 
and  used  as  an  initial  condition  in  the  theory  curves.  The 


theory  that  is  compared  to  the  data  is  also  truncated  at  nmax. 
This  has  only  a  minor  effect  on  the  purity  calculation,  but 
does  change  the  qualitative  result  of  the  entropy,  which  for  a 
limited  basis  set  saturates  much  more  quickly  and  to  a  lower 
value  than  it  otherwise  would.  The  predicted  entropy  and 
purity  curves  (with  no  truncation,  that  is,  as  «max— »°°)  are 
shown  as  the  dotted  lines  in  their  respective  figures. 

C.  Zero-temperature  reservoir 

a.  Decoherence.  Here  we  discuss  the  effect  of  the 
variable-bandwidth  zero-temperature  reservoir  presented  in 
Sec.  Ill  D  4  on  a  superposition  of  Fock  states.  As  before,  we 
prepare  the  superposition  state,  apply  the  reservoir,  and  undo 
the  state  preparation  forming  a  Ramsey  interferometer  to 
measure  the  coherence.  In  this  case,  the  only  state  considered 
was  the  superposition  state  ||)[|0)  +  |2)]/V2. 

We  distinguish  the  two  limiting  cases  discussed  in  Sec. 
Ill  D  4.  In  the  first,  the  effective  decay  rate  yeff  is  much 
faster  than  the  coherent  Rabi  rate  !!lsh .  This  regime  leads  to 
exponential  decay  of  the  coherence  and  is  a  simulation  of  a 
broad  bandwidth  T=  0  reservoir.  In  the  other,  the  effective 
decay  rate  is  slow  compared  to  the  coherent  Rabi  rate.  Here, 
the  system  is  driven  coherently  before  appreciable  decay  sets 
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FIG.  17.  Evolution  of  entropy  for  initial  Fock 
states  in  hot  amplitude  reservoir.  Each  graph 
shows  the  time  evolution  of  the  entropy  of  the 
motional  state  for  an  experiment  starting  in  the 
initial  state  indicated.  The  entropy  is  defined  in 
Eq.  (88),  from  which  the  solid  curves  are  derived. 
The  dotted  lines  are  an  untruncated  theory, 
whereas  the  solid  lines  are  truncated  at  a  maxi¬ 
mum  value  of  n  =  nmax. 


in.  Thus  the  coherent  drive  [on  the  red  side  band  (rsb)]  is 
capable  of  completing  several  Rabi  oscillations  of  the  spin. 
The  evolution  of  fringe  contrast  shows  that  the  coherence 
goes  to  zero  (as  a  7r  pulse  transfers  most  of  the  population 
out  of  1 1 )  1 2 )  into  1 1 )  1 1 )  and  therefore  destroys  all  coherence 
between  |2)  and  |0))  but  then  comes  back  (with  opposite 
sign)  as  a  reservoir  pulse  with  area  greater  than  t t  is  applied. 
During  this  time,  the  coherence  is  decaying  due  to  dissipa¬ 
tion  on  |t).  The  evolution  between  the  two  regimes  is  shown 
in  Fig.  18.  Note  that  the  first  graph  has  almost  no  decay. 
Here  the  red  Doppler  was  absent  ( yeff « IT 0  =  0 )  so  that  the 
residual  decay  is  dominated  by  natural  heating  and  Rabi  fre¬ 
quency  fluctuations.  In  the  subsequent  graphs,  either  the 
strength  of  the  red  Doppler  laser,  or  the  Rabi  frequency  of 
the  rsb  drive,  or  both  were  changed  in  such  a  way  as  to 
increase  the  ratio  of  dissipation  to  coherent  drive,  yeff /Orsb . 

The  experiment  has  an  additional  parameter  that  must  be 
accounted  for.  The  application  of  the  red  Doppler  beam 
causes  a  Stark  shift  of  the  level  |  j~).  In  the  data  of  Fig.  18  this 
effect  was  corrected  for  during  the  experiment  simply  by 
changing  the  detuning  of  the  Raman  beams  to  compensate 
for  the  level  shift;  this  must  be  done  anew  for  each  value  of 
{ I  rsb .  It  is  more  convenient  to  perform  the  experiment  not 


compensating  for  the  Stark  shift  A.  The  theory  for  this  is 
only  slightly  more  complicated  than  that  associated  with  Eq. 
(80),  which  can  be  generalized  to; 


1 

Pm|,0i  =  2eXP(-Tefff/4) 


Oof 


cos- 


5  f lRt 


,  (89) 


~  i  Teff 

A  =  A+-^f,  (90) 

fls=  ^4mg2  +  A2.  (91) 


A  progression  of  increasing  yeff/ns  with  Stark-shift  detun¬ 
ings  is  shown  in  Fig.  19. 

b.  Decay.  For  completeness,  we  also  consider  the  decay 
subject  to  the  zero  temperature  reservoir  of  the  diagonals  of 
the  density  matrix  (with  no  detunings,  A  =  0).  We  measure 
the  time  evolution  of  these  populations  starting  from  an  ini¬ 
tial  state  of  <A(0)  =  1 1)|2).  The  predicted  time  evolution  is 
given  by  Eq.  (12).  The  data  are  shown  in  Fig.  20.  One  will 
immediately  notice  the  negative  “populations.”  This  results 
from  population  left  in  the  |  j )  state  after  the  continuous 


053807-19 


Fringe  Contrast  Fringe  Contrast  Fringe  Contrast  Fringe  Contrast  Fringe  Contrast 


Q.  A.  TURCHETTE  et  al. 


PHYSICAL  REVIEW  A  62  053807 


1 

< 

0.8 

0.6 

0.4 

0.2 

0 

0  10  20  30 

t[ns] 

0. 

0. 

0. 

0. 

0  10  20  30 

t[|0S] 


t[|iS] 


FIG.  18.  Decay  of  Fock  state 
coherences  coupled  to  an  engi¬ 
neered  T=  0  reservoir.  The  points 
in  each  graph  are  the  contrast  of  a 
Ramsey  fringe  that  measures  the 
coherence  between  the  |0)  and  |2) 
motional  states  as  the  time  t  of 
continuous  Raman  cooling  (T=  0 
reservoir)  is  varied.  The  graphs 
have  progressively  larger  ratios  of 
r=yetf/g-  (a)  r=  0.23,  (b)  r=l.l, 
(c)  r  =  2.3,  and  (d)  r=  8.  The 
Stark-induced  detunings  are  cor¬ 
rected  experimentally,  that  is,  the 
position  of  the  RSB  and  its  drive 
laser  are  made  coincident.  The 
theory  curves  are  fits  to  Eq.  (80), 
from  which  r  is  quoted. 


t[|0S] 


t[|iS] 


FIG.  19.  Decay  of  Fock  state 
coherences  in  T=  0  reservoir. 
Same  as  Fig.  18,  except  the  Stark- 
induced  detunings  are  not  cor¬ 
rected  experimentally,  but  are  ac¬ 
counted  for  in  the  theory  curves, 
which  are  fits  to  Eq.  (89).  Front 
the  fits,  yeff/g  =  (0.1,0.7,1.3,2.7, 
2.2)  and  A/g  =  (0,0.5, 0.4,1,2)  for 
graphs  (a)-(e),  respectively.  See 
the  text  for  a  discussion. 
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t[RS] 

FIG.  20.  (Color)  Decay  of  Fock  state  populations  in  a  T=  0  reservoir  starting  from  i/f(0)  =  |2).  The  strength  of  the  red  Doppler  (which 
provides  dissipation)  is  increased  starting  from  the  upper  graph.  Theory  curves  are  from  Eq.  (12)  with  a  modification  to  account  for  residual 
population  left  in  |  )  at  the  end  of  the  continuous  cooling  cycle.  See  the  text  for  a  discussion. 


Raman  cooling.  For  weak  enough  red  Doppler  power,  the 
optical  pumping  does  not  entirely  depump  the  1 1 )  state.  This 
appears  as  a  negative  contribution  to  the  flopping  state  analy¬ 
sis  seen  in  the  first  two  panels  of  Fig.  20  [see  Eq.  (85)]. 

V.  CONCLUSION 

We  have  presented  a  detailed  study  of  the  interaction  of 
several  types  of  quantum  states  with  several  types  of  engi¬ 


neered  reservoirs.  The  results  compare  favorably  with  theo¬ 
retical  predictions,  and  for  several  cases,  we  have  confirmed 
predicted  scalings  for  a  range  of  experimental  conditions. 
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